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ABSTRACT 

A  dimensional  analysis  is  applied  to  the  circulation  around 
Drake  Passage.  It  indicates  presence  of  the  frictional  boundary  layer 
(due  to  horizontal  eddy  viscosity)  with  a  width  of  the  order  of  10  km 
off  the  meridional  coasts,  in  contrast  with  the  frictional  and  inertia 
boundary  layer  with  a  width  of  the  order  of  100  km  in  the  mid-latitude 
ocean.  The  western  boundary  current  off  the  east  coast  of  South  America 
flows  northward  and  southward  respectively  in  the  southern  and  northern 
part  and  has  a  width  of  about  300  km.  The  inertia  terms  tend  to  increase 
the  widths  of  both  the  western  boundary  layer  and  the  current  there.  The 
boundary  current  off  the  west  coast  of  South  America  is  broad  and  ill- 
defined  and  its  intensity  and  direction  is  governed  by  the  wind  stress 
distribution  there  as  well  as  the  bottom  topography  and  the  coastline 
orientation.  The  southern  boundary  layer  off  the  Antarctic  Continent  in 
general  is  much  broader  than  off  the  meridional  coasts  and  has  a  width  of 
the  order  of  600  km  but  in  some  locations  like  Weddel  Sea  and  Ross  Sea 
there  are  jet-like  eastward  currents  due  to  a  leeward  eddy  effect  of  the 
coastline,  forming  counter  currents  there.  The  dimensionless  equations 
of  motion  with  the  inertia  terms  show  that  the  streamlines  are  parallel 
to  the  curves  of  constant  f /H  when  the  terms  with  factors  less  than  10 
are  neglected,  where  f  is  the  Coriolis  parameter  and  H  is  the  depth. 

The  circulation  integral  along  one  of  such  closed  streamlines  indicates 
that  the  transport  divergence  or  upwelling  near  the  bottom  of  the  order 
of  A  x  10-5  (cm/sec)  is  necessary  to  obtain  the  observed  transport 
corresponding  to  the  wind  stress. 
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1.  INTRODUCTION 

It  is  necessary  to  use  approximation  techniques  for  studying  the 
current  in  the  neighborhood  of  the  coasts  of  the  Antarctic  Ocean  and 
particularly  in  the  Drake  Passage,  because  some  terms  neglected  in  (I) 
and  (II)  may  become  important  in  such  areas. 

Equations  of  motion  and  continuity  for  the  volume  transports  in 
the  barotropic  mode  are  given  by 


M  ■  7  ( M  1-r')  -  4  fe  X  r  1  = 


rs  _ i  \ 


.  =± 


i  it  v'V  .m  tAV2  P\  (1) 


’A  =0  (2) 

~> 

where  ii  is  the  transport  vector  with  components  .x  and  liy,  T  is  the 
wind  stress  vector  k  is  the  unit  vector  in  the  vertical  (upward)  direc¬ 
tion  f  is  the  positive  Coriolis  parameter ,  L  and  D  are  the  depth  of  the 
ocean  and  of  the  bottom  Lkman  layer  respectively,  ^  is  the  elevation  of 
the  sea  surface  and  7  and  \Jy  are  divergence  and  Laplacian  in  the  hori 
zontal  level  respectively. 


2.  BOUNDARY  LAYER  DISTANCE 

The  solution  of  equations  of  (1)  and  (2)  can  be  expressed  in  expan¬ 
sions  about  small  parameters  -\) ,  'Moand  1)  =  A/i0  \  ?7>  v;here  H0,  f0 
and  L  are  characteristic  values  of  depth,  Coriolis  parameter  and  hori¬ 
zontal  distance,  respectively.  (II.  Hereafter  the  previous  papers 
"Dynamics  of  A.C.C.*  are  referred  as  (I)  and  (II),  respectively.)  For 
the  A.C.C.  =  2.5  x  1 3* 3  and  )7  =  10 A  The  zeroth  order  term  of  M  in 
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tiiis  expanision  must  satisfy  the  condition  that  the  streamlines  be 
parallel  to  the  contours  of  f/H  =  constant.  T.Jhen  the  coordinate 
system  is  changed  from  the  x-y  system  to  the  ~  ^  system  consisting 
of  the  i  -curves  defined  by  —  ^  =  f /H  and  the  >  /  -curves  orthogonal 

— 7 

to  ^  the  ')  -component  of  the  zeroth  order  term  of  11  should  vanish. 
This  condition  is  not  compatible  with  the  boundary  condition  near  Drake 
Passage,  where  the  streamlines  should  be  parallel  to  the  coastlines. 

There  the  neglected  terms  in  equation  (1)  in  (Il)should  be  taken  into 
account  in  order  to  satisfy  the  boundary  condition. 

A  new  set  of  the  characteristic  values  are  chosen  in  this  section 
to  evaluate  different  terms  of  equations  of  motion  near  Drake  Passage. 

The  characteristic  value  of  the  transport  MQ  (the  zeroth  order  term  in 
the  expansion)  can  be  taken  as  the  mean  value  of  the  zonal  transport  of 
the  A.C.C.  It  equals  2  x  10^  (cm^/sec) ,  because  the  average  of  the  non- 
dimensional  transport  M0  is  about  0.5  (II).  The  characteristic 

horizontal  distance  L  and  surface  elevation  ^  should  be  determined  by 

s  o 

considering  balance  among  different  terms.  In  the  neighborhood  of  the 
coasts  of  South  America  which  are  almost  meridionals  the  vector  ,*i  should 
become  almost  parallel  to  the  coast  or  become  meridional.  Therefore,  in 
the  zonal  component  of  the  equations  of  motion  (1) ,  in  which  the  x-  and 
y-axis  are  taken  eastward  and  northward  respectively,  the  Coriolis  term 
f  iy  has  the  order  of  magnitude  of  fQ  uQ  =  2  x  10^  (cm^/sec2).  The 
Coriolis  terms  in  equation  (A)  should  be  balanced  by  other  terms  as  well 
as  the  pressure  gradient  terms  near  the  meridional  boundaries.  Otherwise, 
XJ  ^  |-| "A.  M  =  q  and  thus  near  the  meridional  boundaries  ily  ^  0 
because  Lix  =0  at  the  boundaries.  The  wind  stress  and  bottom 
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friction  are  smaller  by  two  orders  of  magnitude  than  the  Coriolis  term 
and  the  pressure  gradient  term  respectively  because  T0  =  1  (cm^/sec^)  and 
D/li  =  2.5  x  10~3.  Therefore,  the  advective  and  horizontal  eddy  viscosity 
terms  may  likely  balance  the  Cdriolis  term.  The  former  and  the  latter 
have  orders  of  magnitudes  of  :o2/HL  and  AI10/L^  respectively.  If  the 
advective  terms  balance  the  Coriolis  term,  f  J  =  MQ2/UL,  yielding 
L  =  H0/(fh)  =  5  x  10^  (cm).  This  distance  scale  is  too  small  for  the 

ocean  current.  If  the  horizontal  eddy  viscosity  term  balance  the  Coriolis 
term  f  l‘I0  =  H0/L^s  yielding  L  =  (A/f)'*'^-  10°  (cm).  This  value  may  be 
reasonable  and  is  equivalent  to  the  width  of  the  horizontal  frictional 
boundary  layer.  But  this  width  is  smaller  by  one  order  of  magnitude  than  the 
one  of  the  boundary  layer  along  the  western  coast  of  the  mid-latitude  ocean 
(Stommel,  1963). 

In  the  mid-latitude  ocean,  the  Sverdrup  transport  ly  balances  the 
wind-stress  curl  the  advective  and  Horizontal  eddy  viscosity  terms  of 
vorticity  transport,  where  ^3  is  the  latitudinal  change  of  the  Coriolis 
parameter  Therefore,  the  width  of  the  boundary  layer  cati  be  determined 
from  either  one  of  the  following  relations 


If  we  take  H0  -  2  x  10 '  (cm^/ser.)  and  H  =  2  x  10^  (cm),  1^  becomes  10?  (cm) 

both  from  relations  (3)  and  (4),  in  agreement  with  the  observed  value  of 

width  of  the  western  boundary  layer  (Stommel,  1963).  It  is  to  be  noted 

J- 

that  is  proportional  to  (A/  :  )  ■  or  to  the  one-  third  power  of  the 
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uorizontal  eddy  viscosity  instead  of  A  as  in  the  Antarctic  Ocean. 

Also  the  order  of  magnitude  of  the  advective  terms  represented  by  the 
r.h  s.  of  (4)  is  of  the  same  order  of  magnitdue  as  the  horizontal  eddy 
viscosity  term  in  the  mid-latitude  ocean  when  A  is  taken  as  4000  m.  How¬ 
ever,  the  baroclinic  mode  is  prevalent  in  the  mid-latitude  ocean.  There¬ 
fore  when  the  characteristic .depth  is- taken  as  about  400  m  instead  of 
4000  m  equation  (4)  yields  Lv;  =  224  km.  This  length  seems  to  be  equiv¬ 
alent  to  the  distance  from  the  western  coast  to  the  seaward  boundary  of 
the  western  boundary  current  (the  Gulf  Stream  and  the  Kuroshio).  (If 
the  characteristic  depth  is  taken  as  200  m  or  the  depth  of  the  layer  above 
the  main  therrocline,  the  characteristic  transport  should  be  reduced  to 
half  of  the  assumed  value  that  is  10^  (cm^/sec).  Then  the  value  1^,;  remains 
same. )  The  ratio  of  the  horizontal  eddy  viscosity  terms  to  the  advective 
terms  become  0.09  with  these  values  of  H0  and  Ly.  Therefore  for  the 
baroclinic  mode  the  western  boundary  current  can  be  treated  as  an  inertia 
boundary  lave  ^  /benomenon. 

It  is  now  proved  that  in  the  coastal  boundary  layer  of  the  Antarctic 
Ocean  the  Coriolis  terms,  the  pressure  gradient  and  the  horizontal  eddy 
viscosity  terms  are  balanced  each  other.  The  balance  between  the  first 
two  terms  yields  O  -  0.5  cm.  This  value  is  smaller  by  two  orders  of 
magnitudes  that,  r.hu  v..‘ve  for  the  entire  Antarctic  Ocean,  where  L  is  of 
the  order  of  *0'*  cm.  (  :.I . 


3.  DIMENSIONAL  ANALYSIS  07  .jGUIIDARY  LAYERS 


In  order  to  discuss  the  circulation  near  the  coast  of  the  Antarctic 
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Ocean  the  boundary  layer  techniques  may  be  applied.  The  equations  of 
motion  (1)  can  be  written  in  a  dimensionless  form  as 


-k  K-fn  =  -hv$ 


-> 
T  - 


-f-V  s  ,vfe  +^v  n 

(5) 


where  u  =  i-  /a  is  the  mean  velocity  and  o  is  the  Rossby  number 
(  =M0/fo  L  Kq  )  with  a  numerical  value  of  10"^  corresponding  to  the 
same  characteristic  scale  values  as  in  (II),  where  H0  =  4  x  10^  (cm^/sec), 

L  =  10^  cm  and  other  values  are  the  same  as  in  Section  2.  The  third  term 
of  the  rh.s.  of  equation  (5)  represents  the  bottom  friction  and  can  be 
expressed  by  -  6  f  M  if  the  terms  with  factors  of  or  higher 
orders  of  £  are  neglected  (II).  With  this  bottom  friction  term,  equation 
(5)  can  be  written  as 


a  x  (f  - 


~  ry 

O  /  f 
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~V/£t 
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:>UV 


t  <C 


-e-^M  H H' 


(6) 


wher-  ■?  f s  a  vertical  component  of  the  relative  vorticity  expressed  by 

v7  =  Curl  be  =  3  U  „  /3X  ~  z  1Jr  X  /  '  y 

Z  j 

The  vr ,  i  H  istion  obtained  by  applying  the  operator  Curl  .  to 

equetica  (6)  is 
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where  thr  vn  toner;,  -a  is  defined  by 
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The  dimensionless  parameters  0  and  l''  include  the  characteristic 
values  of  the  vertical  and  horizontal  eddy  viscosity,  respectively.  These 
quantities  are  not  universal  constants  but  may  be  variable.  The  dimension¬ 
less  values  of  the  depth  of  the  bottom  Ekman  layer,  vertical  and  horizontal 
eddy  viscosity  are  designated  by  D,  K  and  A,  respectively.  If  these 
variables  satisfy  the  following  relations 


D  \  H  ^  j  K  II 


1st  (b)  A  H  t'ons?  (C\ 

O  ) 


Co'v 


(9) 


equation  (7)  can  be  written  as 

*  (TAv<P 


V(r/M  H&*V 

=  £  ^  rml  f/H  -  b0Vz(f>)  V  V7'/' 


(10) 


where  the  factors  b  =bf t'A  and  a  =  A/H  are  constant  and  are  assumed  to 

C  0 

be  nearly  equal  to  unity.  The  relative  vorticity  '/  can  be  expressed  with 

/  .  i 

t* 

y  as 

--  V-  (hr1  vft) 


(11) 


'  hen  the  x-y  coordinate  system  is  changed  into  the  £'—^7  system 
in  (II)  equation  (10)  can  be  transformed  into 


as 
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where  and  the  curves  are  orthogonal  to  the  curves. 

The  relative  vorticity  f  and  the  Laplacian  can  ^e  exPtessed 
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respectively  with 
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4 .  WESTERN  BOUNDARY  CURRENTS 

Since  equation  (12)  is  more  suitable  to  aeal  with  the  circulation 
in  the  Antarctic  Ocean  the  coordinates  will  be  used  for  the  boundary 

layer  techniques.  These  coordinates  determined  in  details  near  Drake 
Passage  are  plotted  in  Fig.  1  which  shows  that  the  ^7  curves  for  s? 
larger  than  say  -1.1 ,  are  not  continuous  across  Drake  Passage. 

The  most  likely  site  of  the  western  boundary  current  in  the 
Antarctic  Ocean  is  to  the  east  of  South  America.  The  edge  of  the 
Continental  Shelf  off  the  east  coast  of  South  America  is  the  coastal 
boundary  runs  approximately  normal  to  the  ^7-  curves.  Thus  the 
boundary  is  taken  as  ^  =  Q  .  The  boundary  conditions  for < j  become 

V’  ^  CP'  —  O  <-<  f  -z  o 

Z  '  (15) 

indicating  tnat  both  tangential  and  normal  components  of  the  transport 
should  vanish  along  the  coast.  The  stream  function  p  becomes  finite 
and  has  an  order  of  unity  to  the  east  of  the  boundary  =  0  and  in 
equation  (12)  the  highly  differentiated  terms  have  small  coefficients 
<£  and  }/  .  These  two  conditions  lead  to  the  existence  of  a 
boundary  layer  near  >  =  Q  (Carrier,  1953).  In  this  example  there 
are  three  small  coefficients  instead  of  one  in  conventional  boundary 
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layer  problems.  The^coordinate  is  stretched  by  a  transformation 


}(  _  d  £  —  ^  cl  / 


(16) 


where  s  is  a  new  coordinate  normal  to  the  boundary.  When  only  the 
leading,  terms  are  retained  equation  (12)  becomes 

-  )(~'  j  s /-/  Mg  [(Si1  /d,ss)fi®  A/ 

/  I  /  \  /  .)  ^ 


-  (3  4  /3M)  i  (W 1 3  H ><  v  ^  j  -  fc  P ‘  %  > 2  <S 

«  d>/  •  j  W  o 


■+■  «0V 


/  -  5 


(17) 


where  u  is  the  mean  depth  near  the  boundary  -  In  order  that 

tne  last  term  may  balance  the  first  term,  the  coefficient  of  the  last 
term  should  have  an  order  of  magnitude  of  unity  thus  n  =  1/3.  Then 
the  coefficients  of  the  second  and  the  third  term  of  equation  (17) 
becomes 


3  v 


-Z~h 


/o'  1/3  =-0^6  ;  <z_v  r.rs 


(18) 


Equation  (17)  without  tae  second  term  has  a  solution  of  a  type 


p  =  N  (y)  'J-(A) 


(19) 


where  1W£)  =  Q  ,  =0  and  tj-{c ))  =  1. 


The  function  N  ( )  is 


determined  from  the  inner  solution 


M  -+  t. hr'  ( nip )  A/ /a g  =  o 

'  ( 


(20) 


where  rl  .  is  the  meridional  transport  of  the  first  order  term  in  the 

•  — ^ 

expansion  of  M  about  (II).  If  a  solution  of  equation  (17)  is 
assumed  to  nave  a  form  similar  to  (19) ,  the  dependence  of  /  on  Py 
can  be  parameterized  in  the  non-linear  second  term.  Then  the  ordinary 
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differentia]  equation  for  about  ']  becomes 
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f  l.MiViV/tj  -  /.  f 


(21) 


where  the  prime  indicates  differentiation  with 


and 
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—  f  -> 


2/  xn  hi 'z  =-<■  d-t,  H~l 
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cd 
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1  -  ---  . ? 

The  function  should  satisfy  the  boundary  conditions  that 


(22) 
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In  order  to  obtain  an  approximate  solution  of  (21) ,  the  second  term 
is  linearized  by  substituting  and"’/'  *  with  their  average  values  and 
is  expressed  by 


y  vy  vv )  ~  (V v7 y )  -j- 


/'/ 


(24) 


These  averages  may  be  approximated  by  use  of  a  solution  of  (21)  without 
the  second  term.  For  instance 


/r(  l0^r) 


(25) 


where  the  average  is  taken  for  the  interval  of  (  between  0  and  s  Q 

/ 

and  Q  'z  V  h.  by  use  of  the  substitution  of  (24),  the  solution  of  (21) 

7o  •  y 

satisfying  the  boundary  conditions  (23)  can  be  expressed  by 
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In  tnis  equation  and  m2  are  conjugate  complex  roots  of  the  cubic 
equation 
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where 
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The  numerical  values  of  the  coefficients  in  relations  (28)  are 
determined  for  the  Antarctic  Ocean  between  W  20°  and  W  55°  and  listed  in 
Table  1. 

TABLE  I 


A.L. 

V 

r* 

H6 

•  r3 

a 

b 

-0.70 

34 

- 

5 

1  1  b 

9,1 

0. 8 

1.1 

0.11 

1.24 

0.45 

0,141 

0.363 

-0.75 

39 

6.1 

1.1 

1.4 

0,16 

1.20 

0.60 

0.133 

0.412 

-0.80 

42 

5.2 

1.3 

1.6 

0.19 

1.16 

0.68 

0.131 

0.423 

-0.85 

44 

3.7 

1.5 

1.7 

0.27 

1.10 

0.78 

0.145 

0.443 

0,90 

46 

3.5 

1.7 

1.8 

0.29 

1.04 

0.87 

0.137 

0.487 

-0.95 

49 

Remark  (1)  A.L.  Average  Latitude  of  an  individual  '/  curve 

(2)  Ivi,,  „  is  computed  by  Hpr=  -  T  ?  (^H)  ^ 

5  ~> 

from  equation  (20)  without  the  term  .Lj  r> , 

When  the  roots  m-^  2  are  expressed  with 


.VI 


,9  -  O/f 


(29) 


the  solution  (26)  is  written  as 

-  l  -  e  ^  6b  4  C(  s hs) 


(30) 
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T^e  values  of  a  and  b  are  also  listed  in  Table  1.  The  function  y  which 
is  proportional  to  the  meridional  transport  is  plotted  in  Fig.  2  as  a 
function  of  the  distance  along  -  .  The  values  of  a  and  o  for  this 

curve  is  taken  as  the  average  of  the  values  listed  in  Table  1.  This 
curve  indicates  that  the  western  boundary  current  has  a  width  of  about 
315  km.  The  current  is  broader  and  its  peak  is  less  pronounced  than 
the  one  in  the  middle  latitude  ocean. 

The  direction  of  the  western  boundary  current  depends  on  the  sign 

k  i  - 

of  the  function  \<'  (  ),  However  the  function  /y  (  y  )  can  be  deter¬ 

mined  only  for  its  derivative  with  y  .  The  problem  can  be  solved  by 
considering  the  streamlines  south  and  west  of  the  area  considered. 

.Jhen  streamlines  continue  from  Drake  Passage  to  the  east  of  South  America, 
their  part  corresponding  to  the  western  boundary  current  becomes  north 
ward  When  the  streamlines  of  the  western  boundary  current  are  not 
continuation  from  the  ones  of  Drake  Passage  the  flow  must  be  south¬ 
ward.  This  is  schematically  shown  in  Fig.  3A.  On  the  other  hand  if 
the  stream  function  is  put  to  zero  along  the  coast  of  the  Antarctic 
Continent  as  in  the  previous  studies  (Kamentovich .  1962)  then  the 
western  boundary  current  is  always  northward  to  the  east  of  South  America 
as  shown  in  Fig.  3h.  The  results  obtained  by  Kort  (1959)  and  with 
rotating  model  experiments  (I,  II)  indicate  that  some  streamlines  in 
Drake  Passage  continue  into  the  narrow  western  boundary  current  of  north¬ 
ward  direction  to  the  east  of  South  America  but  also  there  is  a  southward 
boundary  current  due  to  the  streamlines  coming  from  the  east.  Thus  the 
pattern  of  Fig.  3A  is  more  reasonable  than  the  one  of  Fig.  3B.  The 
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streaulines  determined  from  the  solution  (26)  with  numerical  values  of 
a  and  b  listed  in  Table  1  are  plotted  in  Fig.  4. 

In  order  to  study  the  effect  of  the  nonlinear  term  in  equation  (21) 
the  roots  of  equation  (27)  are  determined  by  neglecting  r3  .  The 
real  and  imaginary  parts  of  the  conjugate  roots  defined  by  (29)  are 
shown  in  Table  2. 

TABLE  2 

V  0.70  0.75  0.80  0.85  0.90  0,95 

a  0.258  0.2G8  0.302  0.337  0.3A5 

b  0  382  0.AA3  0.A71  0.536  0.551 

In  Fig.  2  the  function  ij  '  for  the  averaged  values  of  a  and  b  listed  in 

this  table  is  plotted  against  .  Comparison  with  the  curve  for  the 

function  with  the  nonlinear  term  indicates  that  the  nonlinear  effect  in¬ 
creases  the  width  of  the  western  boundary  current. 

The  roots  of  equation  (27)  depend  particularly  on  rp  which  equals 
^ y  I  .  The  values  of  )f  h  are  determined  from  the  distance  of 

excessive  h  curves.  Therefore  the  behavior  of  the  western  boundary 

currents  depends  on  the  depth  distribution  as  well  as  other  physical 
parameters.  If  the  depth  near  the  west  coast  of  South  America  is  constant 
the  values  of  (=  -f/H  )  depends  on  the  latitude  only.  When  the  depth 
'.s  assumed  to  be  5000  m,  the  latitudes  covered  with  lines  ^  =  -0.75 
and  -0.90  approximately  correspond  to  those  covered  by  the  lines  /"  = 

-C.75  and  =  -1.00  for  the  actual  depth  distribution.  The  coefficients 
listed  in  Table  1  are  determined  for  this  fictitious  ocean  of  the  constant 
depth  ana  are  listed  in  Table  3, 
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TAJLE3 


Y 

-0.75 

K.L.  )( f 

40.6 

, '  i  o  c 
> 

-r3 

"rl 

a 

b 

0.80 

7,00 

43.4 

1.20 

1.24 

0.40 

0.143 

0.105 

0.412 

-0.85 

7.42 

46,8 

1.60 

1.55 

0.50 

0.135 

0.057 

0.370 

-0.90 

7.85 

50.6 

1.70 

1.55 

0.50 

0.127 

0.041 

0.360 

-0.95 

8.27 

54.5 

1.30 

1.23 

0.39 

0.121 

0.093 

0.383 

The  function  li'  for  the 

mean  values  of 

a  and  b 

listed  in 

Table  3  is 

plotted 

against 

^  also  in  Fig 

.  2.  it 

is  seen 

that  the  width  of 

the  western 

boundary 

current  decreases  and 

the  decrement  of 

the  amplitude  along 

the 

distance 

from  the  coast 

becomes 

smaller 

for  the 

ocean  of 

the  constant 

depth.  In  other  words,  the  actual  depth  uistribution  intensifies  the 
western  boundary  current  in  the  Antarctic  Ocean. 

5.  EASTERN  BOUNDARY  CURRENT 

The  ^  -/]  coordinates  are  more  complicated  to  the  west  of  South  America 
than  to  the  east  as  shown  in  Fig.  1.  however  in  order  to  discuss  the 
eastern  boundary  current  qualitatively  the  4000  meter  contour  may  be  taken 
as  the  boundary  instead  of  the  seaward  edge  of  the  Continental  Shelf  west 
of  South  America.  This  simplification  may  be  justified  from  the  situation 
that  the  /7  curves  are  bended  sharply  northwards  and  crowded  near  the  4000  m 
contour  but  the  observed  streamlines  (hcrt  1959)  are  southward  off  the 
west  coast  of  South  America.  The  same  boundary  layer  techniques  as  in  the 
previous  section  can  be  applied  to  the  stream  function  with  the  boundary 
conditions  of  (15).  In  the  area  west  of  South  America,  the  solution  of 
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(21)  should  vanish 

for  s  = 

—  .  The 

eastern  boundary 

may  be  represented 

by  s  = 

0  .  The 

constants  corresponding  to 

those  in 

Table  1  are  deter- 

mined  for 

the  area 

west  of 

South  America  and 

listed  in 

Table  4. 

TABLE 

4 

-0.95 

•’Y 

i  ^ 

5 

M 

1  IOC 

"rl 

~r3 

a(=l/2c)  b 

n.i 

0.89 

0.81 

0.090 

2.0 

0.206  0.570 

-1.00 

8.0 

0.93 

0.790 

0.125 

2.1 

0.198  0.604 

-1.05 

7.5 

0.60 

0.49 

0.133 

2.0 

0.123  0.587 

-1.10 

Remark  (1)  The  curves  change  latitudes  sharply  and  thus 
their  average  latitudes  are  meaningless 
(2)i'ci,is  real  root  of  equation  (28) 

Only  the  real  root  of  equation  (2G)  has  a  positive  real  part  which 
makes  an  elementary  solution  vanish  at  the  negative  infinite  s  .  In 
order  that  y  ^ may  vanish  at  s  =  0  ,  the  forcing  function  due  to  the  wind 

stresses  in  equation  (12)  should  be  reinstalled  in  the  r.h.s.  of  equation 

\ 

m  /  , 

(17).  This  term  equals  Si-''  curl  (  ;  i-j Then  the  l.h.s.  of  equation 

V- 

(21)  becomes 
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solution  of  (21)  with  the  term  (31a)  in  the  l.h.s.  becomes 


(32) 


This  solution  satisfies  the  boundary  condition  that  y  =  0  at  s 
Integration  of  (32)  yields 


0  . 


<%  M  i'?)  /  -x $  ■+  (t '  A~i)C-'} 


(33) 


The  function  satisfies  the  condition  that  t /  =  0  at  s  =  0  .  The  condition 


that 


~  £o  M  h/ )  A> 


(34) 


can  be  fulfilled  for  sufficiently  large  \  s\  .  Thus  for  this  range  of 


s  ,  d'  -  |  and  the  function sat;  ■  lies  the  condition  at  s  1. 


The  stream  function  which  is  obtain -.a  by  substituting  (33)  into  (19) 
is  expressed  by 


/  — 


(35) 


Tnis  is  similar  to  the  eastern  boundary  solution  in  the  mid-latitude 
ocean  (Hunk,  1950).  There  is  no  strong  eastern  boundary  current  in  the 
Antarctic  Ocean,  as  seen  from  the  farter  within  the  brackets. 

Since  the  solution  (33)  does  rr<;  lave  an  explicit  term  of  unity  as 
in  (26),  it  is  necessary  to  cbs~k  vrvJ  i  y  of  the  linearizing  process 
applied  to  equation  (21).  It  is  no  a  that  the  exact  solution  of  (21) 
is  expressed  by 
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in  analogy  to  equation  (33),  where  is  dependent  on  ry  and  jJ-^j  j . 

Then  the  second  term  of  equation  (21)  becomes 


I 


'V  — 
I 


<±'"(A  ) 


(37) 


for  a  range  of  s  satisfying  J  sj^.1  ,  when  the  smaller  terms  are 
neglected.  In  a  range  of  negative  s  with  j  s  j  1  satisfying  ^  |  , 

the  nonlinear  term  (37)  can  be  represented  approximately  by  the  third 

order  derivative  .  Therefore,  the  solution  of  equation  (21)  satisfying 

I 

the  condition  of  finiteness  at  0  t  vhas  an  asymptotic  for  exp  (cs  )  for 
this  range  of  s  (Bruijn,  1961). 

The  solution  (35)  indicates  no  jet  1  Ike  eastern  boundary  current. 

If  it  is  assumed  that  the  planetar'7  vcr .  i-:y  (  A  term)  is  balanced  by 
the  bottom  friction  term  only  mat  t.  ‘  tern  boundary,  the  boundary 
layer  equation  (21)  has  only  the  last  ::jo  terms.  The  solution  which 
yields  a  finite  meridional  transport  for-sS^l  is  given  by  the  r.h.s.  of 
equation  (32)  without  the  term  ^  ^ .  In  this  case  the  meridional  transport 
at  the  eastern  coast  s  =  0  is  finite  as  in  the  mid-latitude  ocean  circu¬ 
lation  solved  by  Stommel  (.1948). 

When  the  planetary  vorticity  is  hr,.  ..■ad  by  the  nonlinear  terms  only, 
equation  (21)  with  the  nonhomog°nec»5r.  -  •  can  be  written  as 


This  equation  can  be  transfer:  •' 
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where 


-7 


(AO) 


When  the  r.h.s.  is  neglected  owing  to  its  small  factor  C'Q  ,  equation  (39) 
can  be  integrated  once  with  s  and  becomes 

'f  '  -  ^  (41) 

The  factor  of  the  r.h.s.  of  this  equation  is  an  integration  constant 


chosen  to  satisfy  the  boundary  condition 


The  boundary  layer  solution  of  (41)  is  possible  only  when  .0  .  This 
occurs  if  the  inner  transport  M  is  negative  or  westward.  When  ^  Q  » 
the  solutions  of  the  homogeneous  equation  from  (41)  are  trigonometric 
functions  and  do  not  vanish  for  s  =  f(s>  This  suggests  the  eastward 
current  becomes  unstable  even  near  the  western  or  eastern  boundary  with¬ 
out  frictional  damping  there.  The  boundary  solution  for  is  symme¬ 

trical  about  the  east  and  west  boundaries.  This  solution  corresponds  to 
the  non-viscouse  inertial  current  in  the  mid-latitude  ocean  discussed  by 
Fofonoff  (1954).  The  east-west  asymmetry  of  the  boundary  layer  solution 
will  occur  only  when  the  frictional  force  and  the  small  nonhomogenous 
term  are  taken  into  account.  This  is  similar  to  a  model  of  the  mid¬ 
latitude  circulation  treated  by  Niiler  (1966). 

6.  SOUTHERN  BOUNDARY  LAYER 


The  features  of  the  eastern  and  western  boundary  currents  in  the 
Antarctic  Ocean  are  almost  similar  to  those  in  the  mid-latitude  oceans 
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except  in  characteristic  quantities  such  as  directions  and  intensities. 
The  southern  boundary  formed  by  the  coast  of  the  Antarctic  Continent  is 
unique  in  the  Antarctic  Ocean.  The  boundary  layer  technique  will  be 
applied  again  to  this  boundary. 


The  boundary  layer  equation  in  the  system  can  be  derived  from 
equation  (12)  by  retaining  the  terms  of  the  highest  derivatives  with 
for  each  small  parameter  ■*  ,  ^  or  /  '  .  The  equation  becomes 
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where  the  transformation 


H  >  J  T%  / 9, L 
s4--q 


(42) 


X  J>1  - 


(43) 


is  introduced.  Also  it  is  assumed  that  the  functions  H,  )(  and  }(/} 

§  / 


are  constant  near  the  coast  of  the  Antarctic  Continent.  As  discussed  in 
Report  II,  the  meridional  transport  in  the  Antarctic  Ocean  is  of  the 
order  of  (£  when  it  is  averaged  along  the  whole  latitude  circle.  There¬ 
fore,  if  locally  large  meridional  transports  are  neglected,  the  stream 
function  can  be  expressed  as 
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Then  equation  (42)  becomes 
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If  there  is  a  boundary  layer  near  the  coast  of  the  Antarctic  Continent, 
the  constant  factors  of  the  last  two  terms  in  equation  (45)  should  be  of 


the  same  order  of  magnitude,  because  the  factors  of  other  terms  are  smaller 
by  orders  of  magnitude.  This  leads  to  the  relation  £ -V~"or  i  -  /-?_??  , 
where  £.  -  and  j  .  0.65.  Therefore,  ~ J  ~'V7  r  The  ratios  of 

constant  factors  to  0  of  the  first  to  the  third  terms  are  respectively 

1 ' 


-  ! 


>1 


4-A/c 


(H.J  )( r) 

.-.Tien  the  terms  with  numerical  factors  less  than  10“  1  are  neglected, 
the  last  three  terms  of  equation  (45)  should  be  retained.  Then  the 
solution  of  this  equation  becomes 


d  §c  / V a  - 


for  the  boundary  condition  that 


(46) 


.  /  /  rf  i  -  'A  crt  c 

(«) 

where  1  and  /  is  the  value  of  /  at  s  =  0.  It  is 

°  5  ^  ^ 

assumed  that  the  values  of  j  do  not  change  rapidly  with  increasing  s 

The  boundary  layer  adjacent  to  the  continent  has  a  width  of  ■  /^0>>  = 

628  km  if  the  width  is  defined  as  the  distance  at  which  the  exponential 

terra  of  (46)  reduces  to  exp  (-  '7£  )• 

The  direction  of  the  wind  stress  T_  seems  to  change  from  westward 

to  eastward  within  the  boundary  layer  determined  above  on  the  coast  of 

the  Antarctic  Continent  (Ivanov,  1960).  With  this  value  of  the 

* 

solution  (46)  needs  correction  because  it  is  obtained  under  the  assumption 
that  does  not  vary  rapidly  within  the  boundary  layer.  The  exact  form 

5 

of  the  solution  is  given  by 
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Particularly  when  Tt(.  S  .&*  equation  (48)  becomes 
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The  solution  (49)  indicates  that  there  is  no  ordinary  boundary  layer  in 
which  the  current  changes  rapidly  as  the  increasing  distance  from  the 
coast. 

The  zonal  wind  stress  between  70°  S  and  62.5°  can  be  expressed  by 
the  cosine  function  with  T0  =  -0.6  .  Then  equation  (49)  shows  that  the 
zonal  transport  in  this  latitude  range  is  always  eastward,  notwithstanding 
the  westward  wind  stress.  This  suggests  that  the  easterly  wind  along  the 
coast  of  the  Antarctic  Continent  is  not  a  cause  of  the  counter  current 
near  the  coast.  In  fact,  the  predominance  of  the  last  two  terms  of 
equation  (42)  in  the  boundary  layer  indicates  that  the  lateral  and  bottom 
friction  terms  are  larger  by  one  order  of  magnitudes  than  the  planetary 
vorticity  term,  contrary  to  the  meridional  boundaries.  Therefore,  the 
counter  current  near  the  southern  boundary  is  caused  by  the  frictional 
effect  of  the  coastal  boundary  as  in  case  of  leeward  eddies  behind  an 
obstacle.  This  can  be  confirmed  from  the  observation  that  the  circu¬ 
lation  patterns  obtained  by  the  rotating  tank  experiments  show  the 
westward  current  near  the  Antarctic  coast  in  Weddel  Sea  and  Ross  Sea 
sectors  even  when  the  wind-system  in  the  model  simulates  only  the  t^esterly 
wind  (Report  II).  Mathematically  this  can  be  explained  from  the  boundary 
layer  equation  (45).  Since  }J  a  is  0.2  and  smaller  than  in  the  case  of 
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the  meridional  boundaries,  equation  (45)  should  include  the  second  highest 
derivatives  in  the  bottom  and  lateral  friction  terms.  These  derivatives 
have  factors  containing  the  derivatives  of  X  _  ;t  and  H  which  were 


assumed  to  be  constant  in  derivation  of  (45) .  Then  the  boundary  layer 
equation  becomes 


(50) 


Because  the  coefficients  and  A2  contain  the  derivatives  of  X ^ 
and  E,  they  are  different  in  different  locations  along  the  coastal 
boundary.  The  elementary  solutions  of  the  homogeneous  equation  from  (50) 
are  exponential  functions  of  s  s  the  exponent  m  of  which  satisfies 


I  r  ^ 

a  cubic  equation  similar  to  (27).  When  ^  ~  A,+  77)  A  <0  and 

( A  J  3  A I  --  f\")ZyO+jA^f’  the  solution  of  equation  (50) 


satisfying  the  boundary  conditions  can  be  expressed  with  a  function  simi¬ 
lar  to  (26),  In  this  case,  there  is  a  westward  jet  of  a  narrow  width 
near  the  Antarctic  coast.  This  may  be  a  cause  of  counter  current  near 
the  coast  as  indicated  in  the  observed  surface  current  charts  (Trishnikov, 
1963)  and  in  the  experiments  (Report  II).  It  should  be  noted  that  the 
counter  current  is  local  because  A-^  and  A2  vary  locally. 

7.  CIRCULATION  INTEGRAL  FOR  THE  A.C.C. 

A  mote  general  treatment  of  the  A.C.C.  may  be  carried  out  directly 
by  use  of  equation  (6).  Substitution  of  (8)  into  (6)  yields 
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where 


ch  =  H '' 


7 


(52) 


and 


Z„=  ZH 


I 


J  '  (53) 

The  stream  function  J,  is  expanded  in  a  series  about  small  parameters 
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(r  and  \J  as 
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'ju  ^  i  ' —  'tr' .  *  ’  j:  /  1  t — .  .  ,  ’  (54) 

The  functions  ^  and  should  satisfy  the  following  relations 

(55) 


'N 


/  V #>  =hf'>  (T-fM)-l-/  "(T-  % 4  x  V £  j  (56) 

The  relation  (52)  postulates  that  p  is  a  function  of  £/_  only.  This 
process  suggests  that  the  amplitudes  of  function  ^  cannot  be  determined 
without  taking  into  consideration  the  terms  of  higher  orders  of  £  , 

if  l  and  are  considered  as  independent  of  each  other.  However,  in 
the  Antarctic  Ocean,  and  ^  are  of  the  same  orders  of  the  magnitude 

and  thus  the  terms  C  ~  and  ^cannot  be  treated  as  independent.  This 

‘‘I  O r-' 

indicates  that  the  function  JS  must  be  determined  by  taking  the  non- 

x  / 

linear  advective  terms  into  account.  Near  the  coastal  boundaries,  the 
horizontal  eddy  viscosity  terms  with  parameter  become  important  as 
discussed  in  the  previous  sections. 

In  order  to  obtain  the  relationships  between  the  transport  and  the 
wind  stress,  equation  (51)  is  integrated  along  a  closed  circuit  C.  When 

-t— 

C  is  taken  as  a  curve  on  which  p  is  constant,  O  is  also  constant  on  this 
curve  because  of  (55).  Then  any  single-valued  function  p  the  next  relation 
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is  valid, 

j>c  (57) 

When  the  terms  with  the  second  or  higher  powers  in  and  Jy1  are 

neglected,  the  line  integral  of  equation  (51)  along  C  becomes 

£  ^ i'  -*  V  ^  H'1  lxl  ^  1  =0  ^ 

Equation  (58)  is  similar  to  the  circulation  integral  obtained  by  Niiler 
(1966)  for  the  mid-latitude  oceans.  Equation  (58)  without  the  term  of 
the  parameter  is  similar  to  the  relation  obtained  by  Kamenkovich  (1962) 
for  the  A.C.C..  On  the  other  hand,  in  the  previous  reports  (I,  II),  it 
is  proved  that  the  meridional  transport  of  the  order  of  -  should  be  present 
in  order  to  obtain  a  reasonable  value  of  the  zonal  transport  cf  the  A.C.C. 
for  the  observed  distribution  of  the  zonal  wind  stress.  This  indicates 
that  the  divergence  of  M  does  not  vanish  but  it  has  a  small  value  of  the 
order  of  £  as  expressed  by 

d  /  V  -  £  T)  (59  ) 

-> 

The  transport  vector  II  can  be  divided  into  nondivergent  and  divergent  parts 


as 


where 
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M'l  U  (61) 

Then  the  line  integral  of  the  r.h.s.  of  equation  (6)  on  the  circuit 


C  becomes 
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from  Stokes  theorem,  when  the  terms  with  factors  of  £  or  higher  powers 
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of  ^  and  S  are  neglected.  In  equation  (60),  the  line  integral  along  C  is 
taken  clockwise  eastwards  and  the  l.h.s.  indicates  the  surface  integral 

I  _ / 

with  the  elementary  surface  q  .  The  area  S  is  the  belt  surrounded  by 


the  circuit  C  and  the  Antarctic  Coast  along  which  11  is  assumed  to  vanish. 
When  the  terms  due  to  the  horizontal  eddy  viscosity  are  neglected.,  equation 
(6)  becomes  equivalent  to 
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(61) 


This  is  a  generalized  form  of  Niiler's  relation  (1966)  for  a  circulation 
with  a  small  divergence.  Also  this  equation  was  applied  to  the  A.C.C.  in 
a  special  form  expressed  by  the  coordinate  system  to  determine  the 
meridional  transport  in  Report  II. 

Equation  (61)  indicates  that  the  effect  of  the  meridional  transport 
discussed  in  Report  II  is  equivalent  to  the  divergence.  If  there  is  no 
divergence,  the  magnitude  of  the  transport  of  the  A.C.C.  which  is 
nearly  zonal  should  be  determined  by  equation  (56).  It  was  shown  in  (II) 
that  without  the  divergence,  equation  (56)  yields  the  transport  of  the  A.C.C. 
larger  by  one  order  of  magnitude  than  the  observed  value.  The  divergence 
comes  from  the  thermo-haline  circulation.  This  problem  is  to  be  treated 
in  the  next  report. 

In  order  to  determine  an  average  divergence  in  the  Antarctic  Ocean 
from  equation  (61) ,  the  integration  curve  C  is  taken  as  a  streamline 
corresponding  to  •)  -  -1.27.  This  curve  is  located  at  about  57°  S  as  average 
and  is  almost  along  the  central  line  of  the  transport  belt  which  encircles 
the  Antarctic  Continent  without  being  blocked  by  Drake  Passage  (Kamentovich, 
1962) .  The  numerical  values  of  various  quantities  (in  dimensionless  form) 
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of  equation  (61)  are  as  follows  (II) : 

H  =  0.95,  T  f  »  1.20,  M  =  0.61,  C  -  15.2,  S  =  30.1 
where  C  and  S  are  the  length  of  the  curve  C  and  the  area  between  C  and  the 
Antarctic  Coast,  respectively,  and  M  ^  is  the  average  value  determined  from 
the  observed  transport  obtained  by  Kort  (1959).  Equation  (61)  yields  the 
average  value  of  D  equaling  0.4  or  G  D  ^  10~^.  When  the  dimensional  form 
is  restored,  this  mean  divergence  corresponds  to  the  upwelling  velocity 
of  4  x  10'^  (cm/sec)  near  the  bottom.  The  upwelling  of  this  order  of 
magnitude  is  quite  reasonable. 

8 .  CONCLUDING  REMARKS 

Boundary  layer  techniques  applied  to  the  east  and  west  coasts  of 
South  America  and  to  the  Antarctic  Coast  indicate  the  features  of  the 
boundary  currents  almost  similar  to  those  of  the  mid-latitude  oceans, 
although  characteristic  values  of  the  quantities  representing  the  boundary 
layer  are  widely  different  in  the  two  kinds  of  oceans.  The  cause  of  such 
differences  may  be  partly  due  to  the  predominance  of  a  barotropic  mode 
in  the  Antarctic  Ocean  and  partly  due  to  differences  in  bottom  and 
coastal  topography. 

The  importance  of  the  divergence  of  the  horizontal  transport  to  the 
A.C.C.  is  again  demonstrated  with  the  circulation  integral.  The  positive 
divergence  determined  from  the  observed  data  of  the  wind  stress  and  trans¬ 
port  postulates  the  presence  of  upwelling  near  the  bottom.  Such  process 
may  be  explained  by  considering  the  thermo-haline  circulation  of  the  A.C.C. . 

This  work  was  supported  by  a  grant  NSF  GA-258  from  the  National  Science 


Foundation . 
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EXPLANATION  OF  FIGURES 


Fig.  1  Distribution  of  ^  (  s-f/H)  in  the  neighborhood  of  Drake  Passage, 
Isobaths  of  every  thousand  meters  are  shown  from  Atlas  Antarktiki 
(1)  by  Central  Agency  of  Geodesy  and  Cartography  of  U.S.S.R,, 

Fig,  2  The  values  of  ^  '  against  the  distance  (in  100  km)  from  the 

eastern  coast  of  South  America,  The  curves  I,  II  and  III  are 
for  the  sets  of  a  and  b  in  Table  1,  2  and  3,  respectively. 

Fig#  3  Schematical  distributions  of  streamlines  near  Drake  Passage: 

(A)  Actual  case  (B)  A  case  in  which  the  boundary  value  of 
the  stream  function  is  set  to  zero  along  the  edge  of  the 
Continental  Shelf  off  the  South  American  coast. 

Fig,  1*  Streamlines  determined  by  taking  into  account  the  boundary 
layer  analysis  which  is  summarized  in  Table  1  and  U  for  the 
east  and  west  coast  of  South  America,  respectively. 
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ABSTRACT 


1. 


Dimensionless  forms  of  the  linearized  equations  of 
motion  are  derived  about  the  transports  of  the  A.C.C. 
which  is  assumed  to  be  barotropic.  The  scaling  process 
indicates  that  the  geostrophic  equilibrium  is  predominant 
and  that  the  zonal  wind  stress  is  balanced  with  the  bottom 
friction  and  the  Coriolis  term.  The  horizontal  eddy  vis¬ 
cosity  terms  are  less  than  these  terms  at  least  by  one 
order  of  magnitude.  The  vorticity  relation  postulates 
that  the  streamlines  be  parallel  to  the  isolines  of  f  /  H  , 
where  f  is  the  Coriolis  parameter  and  H  is  the  depth. 
The  magnitude  of  the  quasi-zonal  transport  depends  on  the 
meridional  transport  as  well  as  the  zonal  wind  stress. 

The  average  meridional  transport  determined  from  the 

observed  wind  stress  and  zonal  transport  is  northward 

3  2 

and  about  5x10  cm  /  sec,  almost  equaling  the  value 
estimated  in  the  previous  report.  The  rotating  model 
experiments  indicate  that  the  effect  of  cooling  along  the 
periphery  of  the  Antarctic  continent  generates  eastward 
transport,  in  agreement  with  the  theory.  Also,  the  ex¬ 
periments  show  that  the  main  features  of  the  A.C.C.  are 
maintained  by  the  wind-driven  current. 
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1.  Introduction 


2. 


In  this  paper  bottom  topography,  general  features  of  the  Antarctic 
Circumpolar  Current  (A.C.C.)  in  relation  to  the  effects  of  bottom  topography 
are  discussed.  It  is  assumed  that  the  current  is  barotropic  because  the 
density  stratification  in  the  Antarctic  Ocean  is  small  compared  with  other 
oceans,  as  discussed  by  Kamenkovich  (1962).  The  advective  terms  are 
neglected  because  they  are  smaller  by  at  least  one  to  two  orders  of  magni¬ 
tude  than  other  terms,  as  discussed  in  the  previous  paper  (designated  as  I). 
Although  the  advective  terms  may  become  important  near  the  coast,  their 
effects  may  be  limited  in  the  coastal  boundary  layers.  The  p-plane  is 
used  instead  of  the  spherical  coordinates,  as  used  in  I,  because  the  effects 
of  bottom  topography  are  of  local  nature.  Also,  the  Cartesian  system  is 
easy  to  manipulate.  If  the  meridional  change  of  the  Coriolis  parameter 
is  retained  as  the  p-term,  the  transformation  from  the  Cartesian  to 
spherical  system  is  almost  straightforward. 

2.  Fundamental  Equations 

The  linearized  equations  of  motion  in  the  p-plane  are 


V'  =  -  S*  1 3p/sx  t  t/A^^(i) 


where  the  x-,  y-  and  z-  axes  are  directed  to  the  east,  the  north  and 
the  downward,  respectively  and  f  =  2  £2  cos  0  is  the  Coriolis  parameter 
in  the  southern  co-latitude  0.  The  differentiation  of  f  with  y  yields 
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where  R  is  the  radius  of  the  earth. 


3. 


P  =  -  2  ft  sin  9  /  R, 


When  the  motion  is  assumed  to  be  barotropic,  the  pressure  gradient 


-1 


p  V  p  can  be  replaced  by  the  gradient  of  the  surface  elevation  g  V|<. 
When  equations  (1)  and  (2)  are  integrated  with  z  from  surface  to  bottom 
(the  depth  equaling  H),  the  vertical  eddy  viscosity  terms  can  be  expressed 
by  the  difference  of  the  shearing  stress  at  the  surface  and  at  the  bottom. 
The  surface  shearing  stresses  equal  the  wind  stresses;  while  the  bottom 
shearing  stresses  are  expressed  by 

(K  3U./9Z)  =  -  (Ku)  Jo  :  (Ks  v/az")  —  / D 

hi  iv'  J  y  )"\ 


(3) 


where  D  is  the  thickness  of  the  bottom  boundary  layer  and  and  the  subscript 
H  and  D  refer  to  the  values  at  the  bottom  and  at  the  upper  limit  of  the 
boundary  layer,  respectively  (I).  If  it  is  assumed  that  the  current  becomes 
geostrophic  at  the  depth  H  -  D  (the  upper  limit  of  the  bottom  boundary  layer), 

y  3  3-^/3^  j  VD  —  -f  '%  (4) 

Therefore,  the  bottom  shearing  stresses  become 

(K  a  ul/9z)  ^  =  -g  K  0  (3  ay)  /(D  f)  (5) 

( K  a v-/dZ)H  -  §Kd  05 /«x)/  (0  5-) 

If  the  thickness  of  the  bottom  boundary  layer  is  assumed  to  be  proportional 

,  1/2 

to  the  thickness  of  the  bottom  Ekman  layer,  we  can  put  D  =  (K^  /  f) 

Then  relations  (5)  become 
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which  are  the  same  as  obtained  by  Kamenkovich  (1962). 

The  volume  transport  equations  obtained  by  integrating  of  (1)  and  (2) 
with  respect  to  z  become 


(8) 


are  the 


where  M 


y 


X 


wind  stresses.  In  (I)  the  bottom  shearing  stresses  in  the  x-  and  y-  direc¬ 


tions  are  expressed  with  -kMx  and  -kM  ,  respectively  with  k  »  K  /  DH 


instead  of  the  terms  -g  D  (3£  /  3  y)  and  gD  (dl^/  3  x)  in  Equations  (7)  and 
(8).  The  bottom  stresses  proportional  to  the  transports  were  used  by 
Stommel  (1948)  in  his  model  of  the  circulation  in  a  rectangular  ocean  of  the 
mid-latitude.  The  difference  in  the  result  due  to  these  two  types  of  the 
bottom  friction  terms  is  small,  as  shown  in  the  appendix.  The  integrated 
equation  of  continuity  yields 


0 


(9) 


Equations  (7),  (8)  and  (9)  form  the  system  of  differential  equations  for  un¬ 
known  functions  M  ,  M  and  r  .  Among  the  coefficients  of  the  unknown 
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functions,  Tx,  T^.,  H  and  f  are  considered  as  variable  with  locations,  but 

A  and  D  are  assumed  to  be  constant.  Also,  it  is  assumed  that  D  is  small 

2 

(of  the  order  of  a  few  tens  of  meters)  compared  with  H  and  thus  (D  /  H) 
will  be  neglected  against  unity  in  the  following  derivation. 

Equations  (7)  and  (8)  can  be  transformed  into  non-dimensional  ones  by 
introducing  the  dimensionless  variables  which  are  designated  by  the  primed 
quantities  in  the  following  relations: 

~M  -  Mq  M'  ,  f  =  fQ  f',  H  =  Hq  H',  6  =  60  6',  =  TQ  T' 

(x,  y)  =  L  (  x',  y*  ) 

In  the  Antarctic  Ocean,  the  characteristic  scale  factors  have  the  numerical 
values: 

fQ  =  10“4  (sec-*),  Hq  =  4  x  105  (cm)  L  =  10®  (cm) 


T0  =  1  (cm2/  sec2),  D  =  D0  =  10®(cm)  A  =  AQ  =  10®  (cm^/sec) 


The  characteristic  values  MQ  and  should  be  determined  by  considering 

balance  among  different  terms  of  Equations  (7)  and  (8).  In  (I)  it  is  found 
that  the  zonal  wind  stress  and  the  bottom  shearing  stress  in  Equation  (7) 
balance  each  other.  This  suggests  that  TQ  should  be  taken  equal  to 
g  D  £  Q  /  L,  or  l,o  =  Tq  L  /gD  =  10^  (cm).  In  Equation  (8)  the  geostrophic 
equilibrium  may  be  reached  for  the  meridional  motion  because  the  meri¬ 
dional  wind  stress  T^  almost  vanishes  in  the  Antarctic  Ocean.  Therefore, 
f0  M0  =  gH0  t,  o  /  L  or  M0  =  gH0C0/Lf0  =  T0  Hc/  Df0  =  4  x  106  (cm2/ sec). 
These  values  of  t,  Q  and  MQ  are  the  same  as  obtained  by  Kamenkovich 
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(1962)  through  more  complicated  derivations.  By  use  of  these  character¬ 
istic  values,  the  non-dimensional  form  of  Equations  (7)  and  (8)  becomes 


$  Mu  ~ -  H  /ax  t  £ 


r  c.  9§/ax  fX  V zr\L 


6 


(12) 


where  the  primes  for  variables  are  dropped.  In  these  equations,  the  small 
dimensionless  parameters  <=  and  v  are  defined  by:  e  =  D  /  Hq  = 
2.5  x  10  and  p  =  A  /  fQL  =  10~4.  Equation  (9)  remains  the  same 
in  its  non-dimensional  form. 

The  partial  derivatives  8  l,  / dx  and  8  /  8y  can  be  expressed  with 

wind  stresses  and  transports  by  solving  Equations  (11)  and  (12),  If  terms 
containing  e  ,  €  p  or  higher  orders  of  e  and  v  are  neglected, 

we  have 


3//3X  =  (-  f  n 3  +  £  rxt  U  'v72Mx)H'/  -  si Mx  11  '2(13) 

9^  =  Cihx+£^^V^Ma)l-r'-  H’*  (U) 


Due  to  the  continuity  Equation  (9)  the  transports  Mx  and  M  can  be  ex¬ 
pressed  with  the  stream  function  4  as 


M  x  =  -  d$/3y  ;  M ^ 


9$/3X 


(15) 


Elimination  of  £,  from  (13)  and  (14)  with  substitution  of  (15)  yields 
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H  j(4\  9/ H )  +£H';  {-9  V*'*'  +  (pH /^)  3  ■£>/<?  3 
+  cuyIj/t  +  T*  vH  +  2f  H  (XJH-  V<&)  }  116 

+  ^  1  '-7  4  Vz  $  r  /r  (V^fxVH)j  -  C 

where  J  (  $  ,  f /H)  =  94  /  9x*  0  (f/H)/9y  -  d  $/dyd(f/H)  /  9x  is  a  Jacobian, 
— ^ 

T  is  the  wind  stress  vector  and  the  cross  indicates  the  z -component  of  the 
vectorial  product.  The  quantity  9f/8y  is  dimensionless  and  equals 
-  (sin  0/cos  0Q  )  (  L/R),  where  0  and  0Q  are  co-latitudes  of  the  spot 
considered  and  of  the  reference,  respectively. 

When  the  terms  with  the  order  of  e  or  smaller  values  are  neglected, 
the  stream  function  must  satisfy 

T  (<?  ,  9/ H )  —0  in) 


or  4.  is  a  function  of  f /H  only.  Therefore,  when  H  is  constant,  4  is 
a  function  of  the  latitude  only,  yielding  a  purely  zonal  current.  This  is  a 
Taylor-Proudman  principle  (Proudman,  1953).  The  condition  (17)  is  also 
applied  to  the  A.C.C.  by  Kamenkovich  (1962).  However,  this  does  not  conflict 
with  the  model  proposed  in  (I)  where  it  is  proven  that  the  meridional  compo- 
nent  of  transport  plays  an  important  role.  In  fact,  the  meridional  transport 
in  (I)  is  of  the  order  of  e  and  thus  Equation  (17)  is  not  violated  in  the 
model  of  (I),  although  the  meridional  transport  in  (I)  represents  the  zonally- 
averaged  values*  (When  the  meridional  transport  is  finite  even  though  very 
small,  the  resulting  streamlines  in  a  circular  ocean  like  the  Antarctic  Ocean 


become  spiral). 
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3.  Solution  of  the  basic  equations 


The  solution  of  Equation  (16)  must  satisfy  relation  (17),  which  postulates 
that  4  is  a  function  of  f/H  only.  However,  this  equation  does  not  speci¬ 
fy  the  function  itself.  In  order  to  determine  the  function,  the  sum  of  the 
terms  with  a  factor  e  in  Equation  (16)  should  be  put  to  zero.  The  terms 
with  a  factor  v  may  modify  the  behaviors  of  the  solution  near  the  boundaries, 
It  is  more  convenient  to  use  the  orthogonal  curvilinear  coordinates 


(  4, 


),  in  which  r\  is  defined  by 


7  = 


f/H 


(18) 


and  ^  *s  orthogonal  to  r\  .  Equation  (16)  can  be  written  in  a 
vectorial  form  as 


l 


diir  (  f M  /H )  +  £  1  Ourlz  (T/H)  - 

+  V  Cma J,  (H_/  VZM)=0 


(19) 


where  M  is  the  vector  of  transport  and  curl  is  the  z-  component 


of  the  rotation  of  a  vector.  When  the  coordinate  system  is  changed  from 
the  (x,  y)  coordinates  to  the  (  r\  )  coordinates,  the  divergence  and 


curl  can  be  expressed  with  the  (£,  -n  )  system  as 
z  * 

chv  V/  -  ^ 


l/7  Vj 


(20) 


QMA 


V  =  Ky) 
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d : 


(y  k  )  -  -2-  ( 

\  *i )  3^  \ 


v?  np 


(21) 


where  V,  and  V  are  the  6-  and  -n  -  component  of  the  vector  V, 

4  ^ 


respectively.  The  scale  coefficients  )({ 


){  are  defined  by 
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K 


9 


(22) 


where  V  indicates  the  gradients  of  £  and  -q  as  functions  of  x  and  y. 
By  use  of  the  relation  (20),  Equation  (9)  becomes 


3  (Mr  )  4) 


£ 


c) 


(  !  V/  A?  ^  <T  .)  /  a'  ^7  -  0 


(23) 


The  stream  function  e  in  the  (  | ,  ri  )  coordinate  system  is  defined  by 


M-  =  -  a 


9  7  ; 


a ' 


-/ 


3  3>/S> 


(24) 


By  use  of  (18)  and  (24),  we  have 


dor  (-fH/H)-  (his  (-'/  lvl)~  7v  Ks)  3P/c]5 


(25) 


Equation  (18)  can  be  transformed  into 


O  ic)  ^ 
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-  ?  \JL-\nH  Ifh 

(  9  ?  \  hi  X^ 


3  /  J±Ai 


9^ 


HXf  3', ) 


9  7  \  H 

?  in**  $£. 


-  u 
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(a  v  U 
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(26) 
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by  use  of  relations  (21)  and  (25).  In  this  equation  the  Laplacian  V  can 


be  expressed  by 
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X r  O 


-i  7  a  ( Khj  ) 

l  ^5  Uf 


3  ^  ) 


(27) 


The  solution  of  (26)  can  be  expressed  in  terms  of  expansion  about  4^  and 
y  .  Since  y  is  smaller  by  at  least  one  order  of  magnitude  than  6  »  the 

terms  with  v  may  be  dropped  except  close  to  the  boundaries,  where  the 
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gradients  of  4  perpendicular  to  the  boundaries  become  large. 
4  is  expressed  by 


Then 


10. 


<3?  =  +  8  T  + 


(28) 


Substituting  (28)  into  Equation  (26)  and  equating  the  coefficients  of  the 
equal  power  of  e  to  zero,  we  have  at  once 


(29) 


This  is  another  form  of  relation  (17).  The  same  result  was  obtained  by 
Kamenkovich  (1962),  although  in  his  coordinate  system  £  and  rf  corres¬ 
pond  to  respectively  ^  and  £  in  the  present  system.  The  coefficients 
of  the  first  power  of  e  yields  the  following  equation: 

-  S  9  (T,  Xnl-T)  d(7>  X(  II -j l  3  {  n  K<r  d£c\0 

t  £  *?  c)  j  (  j-j  K/^  3  /(?  / 


(30) 


This  equation  is  also  derived  by  Kamenkovich  (1962),  except  the  first  term 
which  he  neglected.  The  first  term  represents  the  effect  of  the  meridional 
transport,  as  discussed  in  the  previous  report  (I).  This  term  substan¬ 
tially  modifies  the  magnitude  of  the  zonal  transport  due  to  the  zonal  wind 
stresses  only.  Equations(28),  (29)  and  (30)  indicate  that  the  zonal  compo¬ 
nent  of  the  transport  in  the  Antarctic  Ocean  is  of  the  order  of  e  but  it 
plays  an  important  role  to  change  the  wind-driven  zonal  transport.  The 
function  4  cannot  be  determined  from  Equation  (26)  without  considering 
the  terms  with  the  factor  p  and  with  higher  orders  of  e  This  indi¬ 
cates  that  the  function  4  depends  on  the  processes  near  the  boundaries. 
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Also,  the  convection  process  at  the  periphery  of  the  Antarctic  Continent 

may  generate  the  meridional  transport  as  discussed  in  (I), 

It  is  necessary  to  scrutinize  Kamenkovich's  (1962)  results,  because  he 

14  3/ 

obtained  the  total  eastward  transport  of  1.4  x  10  cm  /  sec  without  taking 
into  account  the  meridional  transport.  This  value  is  in  good  agreement 
with  the  observed  value  by  Kort  (1959)  of  1.  3  x  10^  cm^/sec,  but  the 
result  of  the  present  study  (I)  indicates  that  the  eastward  transport  becomes 
higher  without  the  meridional  transport.  A  close  check  of  his  calculation 
reveals  two  reasons  for  his  good  agreement  with  the  observed  data.  One 
reason  is  that  he  chose  the  zoneof  the  A.  C.C.  between  the  isolines  of  4^ 
and  4,5*  the  zonally-averaged  latitudes  of  which  are  respectively  at  59.8°S 
and  53.  3°S,  Therefore,  the  width  of  the  A,  C.C.  in  his  calculation  is  about 
6.  5°  latitudes.  The  other  reason  is  that  his  values  of  F'1'1  in  Table  1 
seem  to  be  inexplicably  low.  This  function  corresponds  to  the  zonal 

average  of  cos  6  *T  -  /H  in  the  notation  of  this  paper,  where  T  and 

YL  a, 

H  are  non-dimensional  zonal  wind  stress  and  depth,  respectively  and  6 
is  the  angle  between  the  n  curve  (  4  curve  in  Kamenkovich's  notation) 

and  the  latitude  circle.  The  result  of  recalculation  of  F  O)  by  use  of 

2 

the  wind  stress  data  of  Ivanov  (I960)  and  bottom  topography  and  4"curves 
of  Kamenkovich  (1962)  shows  that  its  values  at  the  southern  and  northern 
boundaries  of  the  A.  C.  C.  are,  respectively  0.  86  and  - 1.  3  compared  with  - 
-0.46  and  -0,71  of  Kamenkovich.  The  wind  stress  T^  between  60°  S 
and  50°  S  is  above  1.  3  (crrm/ sem)  and  the  depth  is  between  3000  m  and 
4500m,  yielding  H  nearly  equaling  unity.  Even  when  it  is  assumed  that 
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the  average  angle  6  is  15°  and  the  average  depth  is  about  1000  m,  the 
values  of  cos  6*T~  /H  exceed  one.  Therefore,  our  estimate  of  F  ) 

a  ? 

seems  to  be  correct.  , 

Equation  (30)  can  be  used  to  determine  the  zonal  transport.  When 
this  equation  is  average  with  respect  to  4  it  becomes 
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(31) 


where  the  bar  indicates  the  mean  with  respect  to  4  anci  ^  is  the 
meridional  component  of  the  transport  corresponding  to  the  stream  func¬ 
tion  4  .  .  A  boundary  condition  is  that  4  =  o  at  the  southern 

1  o 

boundary  ^  =  r\  . 

s 

Equations  (13)  and  (14)  can  be  written  in  the  vector  notation 

QyoJ  5  =  J/f  rj  ^  e  A  +i^  y27i)fi~>  1 

-  £  f  M  H"2 

where  k  is  a  unit  vector  in  the  vertical  direction. 

When  Equation  (32)  is  expressed  with  the  coordinate  system  (  4»  "H  )»  the 
4-  and  tj  -  component  of  the  equation  become,  respectively 


where  the  Laplacian 


can  be  expressed  by  (27)  and  relation  (18)  is 
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substituted.  When  the  transport  vector  M  is  expanded  in  a  power  series 
with  respect  to  €  ,  v  _ ^ 

M  =  Mq  +  £  Mj  + 

as  in  (28),  the  component  MQ  vanishes  as  seen  in  Equation  (29).  Then 
the  average  of  (33)  with  respect  to  £  becomes 


=  o 

(35) 


If  v  It  is  neglected  against  one,  we  have 


(36) 


Equation  (36)  is  essentially  the  same  as  the  integrated  form  of  Equation  (31) 


with  respect  to  ^  .  When  M  ,  ^  vanishes,  Equation  (36)  yields 
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Therefore,  an  integration  constant  of  Equation  (31)  without  M  ( ^  term 
can  be  determined  to  be  zero,  as  derived  by  Kamenkovich  (1962)  in  a  more 


circuitous  way.  If  it  is  assumed  that 
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Equation  (36)  yields 
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Therefore,  the  zonal  transport  between  A?  —  7?  and  ^  can  be  given  by 
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Further,  it  is  assumed  that  M,  ,  H  -  M  /  .  H  and  that  M  ,  ^  is 

'  t  £  c 


constant  between  f?  ^  and  ^  Then,  when  the  zonal  transport  between 
the  two  streamlines  t)  =  /£  ^  and  ^  ^  is  known  from  observed 


data,  for  instance  by  the  dynamic  calculation,  the  average  value  of  M  t  r 


? 


can  be  determined  from  Equation  (40)  because  the  term  T-.  can  be  deter- 
mined  from  the  zonal  wind  stress.  For  actual  computation,  Kamenkovich' s 
(1962)  streamlines  from  O  =  h  s~  ^  \  to  ^  -  O  n  =  ^  5  anc* 

their  curvature  KA^  are  used(in  his  notation  and  X  ~  »  respectively). 
The  wind  stress  term  T^  is  recalculated,  because  his  calculation  of  the 
wind  stress  term  seems  to  be  erroneous.  The  numerical  values  of  the 
variables  necessary  for  the  present  calculation  are  listed  below: 
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-1.39 

-1. 33 

-1.27 

-1.21 

-1.  17 

2.74 

3.  02 

3.  16 

3.29 

3.  43 

~r ' 

1  5 

-0.86 

-1.05 

-1. 20 

-  1.27 

-  1.  34 

H 

0.90 

0.92 

0.95 

0.97 

1.00 

By  use 

of  these  values, 

the  eastward  transport  due  to 

the  wind 

stress  be- 

tween 

7  -  7  s 

and 

7 =  7  n 

becomes  0.  625. 

The  zonal  transport 
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1 . 3  x  10 ^  (cm^/ sec)  obtained  by  Kort  (1959)  becomes  0.  325  in  the  non- 

dimensional  scale.  Since  the  integral  j  ;/  H  X^  cl  ^  becomes  0.64,  the 

average  value  of  M  equals  0.47  and  its  direction  is  northward.  The 

dimensional  value  of  this  meridional  transport  p  M  ,  x  4  x  10^ 

3  3  3 

(cm  /  sec)and  equals  5  x  lO^  cm^/sec,  This  val  ue  is  nearly  equal  to  the 
estimated  northward  transport  7  x  10"5  (cmJ/sec)  in  Part  I. 

4.  Comparison  with  experimental  results 

In  order  to  study  the  effects  of  cooling  along  the  periphery  of  the 
Antarctic  Continent,  a  series  of  experiments  were  made  with  the  scale 
model  on  a  rotating  table  (Ichiye,  1966).  A  polyethylene  tube  was  filled 
with  ice  cubes  and  installed  along  the  periphery  of  the  model  Antarctic 
Continent.  The  table  was  rotated  with  speeds  of  4,  6  and  8  r.p.m. 
Temperature  was  measured  with  mercury  thermometers  every  few  minutes 
at  several  points  at  the  rim  of  the  circular  tank  and  of  the  edge  of  the 
continent.  The  temperature  differences  between  the  continental  edge  and 
the  rim  of  the  tank  were  usually  about  1*5°  C  for  these  experiments  and 
could  be  maintained  as  such  for  more  than  thirty  minutes.  The  temperature 
differences  at  different  points  along  the  inner  edge  or  along  the  rim  were 
less  than  0.2°  C.  The  16  mm  motion  pictures  were  taken  at  the  rate  of 
four  frames  per  rotation  with  a  camera  rotating  synchronously  with  the 
table.  Small  paper  patches  were  strewed  at  the  surface  to  determine 
surface  circulation.  Ink  spots  were  introduced  at  several  places  to  obtain 


the  bottom  circulation. 


'  i  -  ■  V 

r 

•  '  .  :  .  .  .  _  •  i  - 

'i  - ■ 

.1  >'  b,  .  :  •  •  ■  '  f  .  ''  '  *  ’ 

•  •  ,  I  '  .  ■  J  -  :  Y  C- 1  .  ■  'a*  ..  . 

.  •  .  .  ■  rn  ^  '  ■  ::  \ 

.  1  '•  " .  ■  '  1  .  ■  ■  «  .  •  •  ■  •  /  • 


•  •  A  a 


•  ;  ;-’i  .  ;  •  ■  <  _  • .  '  . 

' 

o  ■  •  ...  ■  •  ' 

v  i:'  ,<i-  /  y  .  •.  .  £  v 

1  ;iv>  '  .  •  v  ’  "  •  is  ■  •  -  ■  •  ■ 

••  *  '  ii.<  •  .  ::  .. . 

»r  <4  f.  i-  •  S'  •  s  A .  i 

t  o  . 

‘  it  u  i  :*x i  :  r.  r.im  0  I  >t 

• :  .  '•  .»  .  '  •  •  T  • 

•  :  •  • .  .  ' .  ■  •  .  .  . 

.  ..  .  . 


16. 


In  Figure  1  one  of  the  circulation  patterns  determined  by  such  films 
is  shown.  The  rotation  rate  was  6  r.p.m.  The  surface  and  bottom  currents 
are  indicated  with  solid  and  dashed  lines,  respectively.  One  of  the  most 
striking  features  of  the  circulation  is  a  general  eastward  flow,  contrary  to 
the  actual  A.C.C.  Particularly  a  strong  current  of  an  easterly  direction 
flows  through  the  Drake  Passage.  This  is  always  the  case  regardless  of 
the  rotation  speed.  Meridional  components  of  the  surface  currents  are 
generally  southward  and  those  of  bottom  currents  are  the  reverse  except 
near  the  Drake  Passage  where  both  are  almost  similar  and  northward  west 
of  it.  The  difference  of  the  meridional  components  between  the  surface 
and  bottom  is  less  for  8  r.p.m.  and  more  for  4  r.p.m.  than  for  6  r.p.m. 

In  the  Weddell  Sea  there  is  a  clockwise  eddy  in  this  experiment.  However, 
experiments  for  8  r.p.m.  showed  a  counter-clockwise  eddy  there.  Also, 
some  details  of  the  circulation  are  different  for  different  sets  of  parameters 
like  the  rotation  speed  or  the  temperature  difference  between  the  continental 
edge  and  the  rim.  Notwithstanding  all  such  details,  the  general  features  of 
Fig.  1  are  closely  in  agreement  with  a  theoretical  prediction  in  Section  III. 
The  cooling  along  the  continental  edge  causes  sinking  there.  The  water 
which  sinks  spreads  north  along  the  bottom  and  then  flows  eastward. 

For  comparison,  one  of  the  patterns  from  the  experiments  for  the 
wind-driven  circulation  is  shown  in  Fig.  2.  The  wind  was  blown  over  the 
model  with  six  wind-blowers  installed  at  the  edge  of  the  table.  The  rotation 
speed  is  again  6  r.p.m.  This  figure  indicates  that  the  circulation  is  similar 
to  the  one  determined  by  Kort  (1959)  and  Treshnikov  (1963)  in  several  points. 
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17. 

Particularly,  the  streamlines  follow  the  bottom  topography  with  a  tendency 
that  they  deviate  southward  over  a  trench  and  northward  over  a  ridge.  Also, 
the  current  along  the  periphery  of  the  continent  is  generally  eastward  and 
there  is  a  clockwise  eddy  in  the  Weddell  Sea.  This  suggests  that  the  A.C.C. 
is  maintained  for  the  most  part  by  the  wind,  as  discussed  in  Part  III,  On 
the  other  hand,  the  effect  of  the  thermohaline  circulation  is  to  counteract 
the  wind  stress  as  demonstrated  in  Fig.  1  and  to  change  the  circulation 
quantitatively . 


•••  •  .  .  :i  '  '  •  '  •  ;  *v  .  v  . 

•  rt 

i  , 

■ 


EXPLANATION  OF  FIGURES 
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Figure  1 


Figure  2 


Streamlines  for  the  cooling  and  heating  experiment 
(R-tation  Rate:  6r.p.m.  Solid  lines:  surface  currents 
determined  with  paper  patches  Dashed  lines:  bot¬ 
tom  currents  determined  with  ink  spAs 

Surface  streamlines  for  the  wind-driven  circulation 
experiment  (Rotation  rate:  6  r.p.m.) 
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ABSTRACT 


Equations  of  motion  and  continuity  of  the  volume  transport  of  the 
Antarctic  Circumpolar  Current  (A.C.C.)  are  derived  for  the  barotropic 
mode,  including  the  advective  terms,  bottom  friction  and  horizontal 
stress  terms.  In  the  zonal  component  of  equations  of  motion  averaged 
along  the  entire  latitude  circles,  the  Coriolis  term,  the  zonal  wind 
stress  and  the  bottom  friction  are  of  the  same  order  of  magnitude,  while 
the  advective  terms,  the  pressure  gradient  and  the  horizontal  eddy  vis¬ 
cosity  terms  are  less  by  one  to  two  orders  of  magnitudes  than  these 
terms.  The  meridional  transport  changes  the  zonal  transport  through 
the  Coriolis  effect  in  such  balance  of  forces  in  the  zonal  direction. 

The  northward  transport  of  about  7  x  leP  (cm^/sec)  which  was  postulated 
by  Munk  (1966)  from  conservation  of  salinity,  carbon- lU,  oxygen  and  heat 
will  decrease  the  eastward  transport  due  to  the  zonal  wind  stress  only 
almost  to  half.  Without  the  effect  of  the  meridional  transport,  the 
zonal  transport  due  to  the  zonal  wind  stress  becomes  almost  double  the 
observed  value.  The  horizontal  eddy  viscosity  terms  becomes  important 
only  near  the  north  and  south  boundaries  of  the  A.C.C.  The  meridional 
transport  seems  to  be  generated  by  sinking  of  the  dense  water  along  the 
periphery  of  the  Antarctic  Continent  as  demonstrated  with  a  rotating 
table  experiment.  The  effect  of  Drake  Passage  also  is  to  decrease  the 
zonal  transport  by  increasing  the  average  northward  transport  as  ver¬ 
ified  in  the  Antarctic  model  experiment  with  a  rotating  table. 
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1.  INTRODUCTION 

The  Antarctic  Circumpolar  Current  is  very  unique  compared  with 
other  current  systems  in  the  ocean.  First,  it  is  not  impeded  by  land 
barriers.  Second,  its  density  stratification  is  less  pronounced  at  all 
the  depths  than  in  the  mid  and  low-latitude  oceans  and  particularly  it 
lacks  a  strong  thermocline  as  found  in  other  oceans.  Third,  the  perma¬ 
nent  wind  field  is  almost  zonal  and  westerly  except  over  the  Antarctic 
Continent# 

There  are  a  number  of  theoretical  works  on  dynamics  of  the  Antarctic 
Circumpolar  Current  (A,C,C,),  These  include  Munk  and  Palm&n  (195>l), 
Hidaka  and  Tsuchiya  (1953),  Fofonoff  (19!?1|.),  Takano  (195&),  Stommel 
(1957),  Wyrtki  (i960),  Kamenkovich  (1962)  and  Barcilon  (1966),  Only 
Fofonoff  and  Barcilon  discussed  the  thermohaline  and  convective  circu¬ 
lation  of  the  current.  Other  works  treated  the  current  as  a  wind- 
driven  circulation  and  are  critically  reviewed  by  McKee  (1966)# 

Earlier  models  of  Munk  and  Palmen,  Hidaka  and  Tsuchiya,  and  Takano, 
are  based  on  the  circular  symetrical  geometry  of  the  Antarctic  Ocean 
with  constant  depth  and  they  obtained  the  zonal  water  transport  with 
larger  by  orders  of  magnitudes  than  the  observed  value  except,  when 
extraordinary  large  values  of  horizontal  eddy  viscosity  are  used, 

Stommel  (1957)  suggested  that  such  discrepancy  may  be  remedied  when 
assumption  of  zonal  uniformity  of  the  current  is  abandoned  as  the  cur¬ 
rent  is  really  not  zonal,  Kamenkovich  (1962)  considered  the  vertical 
eddy  viscosity  terms  only  as  frictional  forces.  Taking  into  account  of 
depth  variation,  Vie  obtained  the  transport  agreeing  with  the  observed 
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value.  Barcilon  (19 66)  suggested  that  the  land  drainage  from  the 
Antarctic  Continent  flows  in  the  reverse  direction  to  the  A.C.C. 
and  thus  reduces  the  transport  of  the  wind-driven  eastward  flow. 

All  the  mathematical  models  are  based  on  simple  bottom  and  coastal 
topography,  because  otherwise  the  equations  become  untangible.  There 
are  two  approaches  to  further  study  of  the  dynamics  of  the  current. 

One  is  analogue  techniques  including  hydrodynamic  experiments  with  a 
scale  model,  electrolytic  tank  and  membrane  methods.  The  other  is 
numerical  analysis  of  the  basic  equations.  The  present  work  is  a 
preliminary  for  interpretation  of  the  experiments  with  a  rotating 
table  which  have  been  carried  out  by  the  author  since  196£. 

2,  BASIC  EQUATIONS 

Since  the  A.C.C,  envelops  the  entire  latitude  circles,  the 
spherical  coordinates  seem  to  be  more  appropriate  than  the  Cartesian 
coordinates.  However,  some  local  features  such  as  the  effects  of 
bottom  topography  or  coastlines  may  be  explained  in  less  complicated 
ways  by  use  of  the  Cartesian  coordinates  as  will  be  discussed  in  the 
second  report. 

The  horizontal  components  of  equations  of  motion  of  the  stationary 
ocean  currents  in  the  surface  spherical  coordinates  are  as  follows: 
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where  £  (j  ^and  the  associated  velocities  w,  v,  u,  correspond  to  the 
directions  along  the  vertical  (radius)  co-latitude  (7F/^„  minus  lati¬ 
tude)  and  longitude  respectively,  R  is  the  radius  of  the  earth,  and  K 
and  A  are  vertical  and  horizontal  eddy  viscosity  respectively.  The  co¬ 
latitude  and  the  longitude  are  measured  respectively  from  south  to  north 
and  from  west  to  east.  The  operator  is  the  horizontal  Laplacian  de¬ 
fined  by 
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(3) 


The  vertical  component  of  the  equations  of  motion  is  simply  given  by  a 
hydrostatic  equation 
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The  equations  of  continuity  and  conservation  of  mass  are  respectively 
given  by 
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3.  ESTIMATION  OF  VARIOUS  TERMS 

Kamenkovich  (1962) estimated  orders  of  magnitude  of  various  terms 
in  the  linearized  equations  of  motion  and  concluded  that  only  the 
Coriolis^ pressure  gradient  and  vertical  eddy  viscosity  terms  are  im¬ 
portant,  It  is  necessary  to  revise  his  estimate  if  we  are  concerned 
with  details  of  the  current  pattern  rather  than  the  gross  feature  of 
the  current. 

When  the  characteristic  values  of  horizontal  and  vertical  distances 
and  horizontal  and  vertical  velocities  are  taken  as  L,  H,  U  and  W  respec¬ 
tively,  the  ratios  of  horizontal  and  vertical  convective  terms  and  hori¬ 
zontal  and  vertical  eddy  viscosity  terms  to  the  Coriolis  term  in  equation 
(1)  and  (2)  are  respectively  expressed  by 


where  y  is  the  characteristic  value  of  the  Coriolis  parameter.  The 
equation  of  continuity 


For  the  Antarctic  Circumpolar  Current,  z  10"*^  (sec“  !),  0  =  10 
(cm/sec),  L  »  10®  (cm)  and  H  =  10 5  ( cm).  Therefore^  W  ^  lCr2  (cm/sec)* 


The  horizontal  and  vertical  Rossby  numbers  R^  and  Rv  are  of  the  orders 
of  10“3#  Thus  the  non-linear  terns  may  be  neglected  in  treating  the 
gross  features  of  the  current.  The  horizontal  and  vertical  eddy  vis¬ 
cosity  A  and  K  can  be  estimated  from  the  empirical  formula  A  =  k^_ 
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and  K  s  k2  respectively,  where  k^  and  k2  are  dimensional  constants 

with  values  ranging  0,01  to  0,02  (Olson  and  Ichiye,  19!?9)  for  L  and  H 
being  expressed  in  cm.  Then  the  horizontal  and  vertical  Ekman  numbers 

* 

and  E^.  are  expressed  by  k^  L “2/3  ^  f-1  K-2/3  respectively, 

yielding  respectively  the  numerical  values  of  10**3  and  lG*-1  for  k^  s 
kg  =  0,02,  These  estimates  indicate  that  the  non-linear  terms  and  the 
horizontal  eddy  viscosity  terms  are  of  the  same  orders  of  magnitudes  and 
may  be  neglected  compared  with  the  Coriolis  terms  and  that  the  vertical 
viscosity  term  may  be  not  negligible  as  concluded  by  Kamenkovich  (1962), 
However,  in  this  estimation  the  li/3-power  law  yields  the  horizontal 
and  vertical  eddy  viscosity  of  atout  10^  and  10^  (cm^/sec)  respectively,, 
compared  with  the  generally  accepted  values  of  10^  and  10^  (cm2/sec)  in 
the  area.  If  the  latter  set  of  eddy  viscosity  is  used,  the  Ekman  numbers 
Eh  and  Ev  are  10"^  and  10-3  respectively.  Therefore,  the  horizontal  and 
vertical  eddy  viscosity  terms  are  respectively  smaller  by  one  order  of 
magnitude  than  and  of  the  same  order  of  magnitude  as  the  advective  terms. 
On  the  other  hand,  the  eddy  viscosity  terms  become  only  important  at  the 
boundary  layer,  the  thickness  of  which  in  turn  depends  on  turbulence  in 
the  current  or  eddy  viscosity.  This  indicates  that  there  is  no  accurate 
method  of  estimating  eddy  viscosity  terms.  In  order  to  obtain  rough 
estimates  of  the  eddy  viscosity  terms  x-rithin  the  boundaxy  layer,  some 
observed  data  must  be  used,  although  these  data  are  variable  at  different 
locations  and  situations.  The  thickness  of  the  surface  and  bottom  bound¬ 
ary  layer  of  the  ocean  ranges  from  10  m  to  100  m,  yielding  the  vertical 
Ekman  number  of  1  to  10"^  for  the  vertical  eddy  viscosity  of  10^  (cm2/sec). 


' 
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The  distance  of  the  coastal  boundary  layer  is  less  certain  than  the  thick¬ 
ness  of  the  surface  and  bottom  boundary  layer  but  it  has  an  order  of  magni¬ 
tudes  of  100  km.  The  horizontal  Ekman  number  E^  becomes  10“ 2  for  this 
distance  corresponding  to  A  s  10®  (cm2/sec).  Then  E^  is  of  the  same  order 
of  magnitudes  as  the  Rossby  number  for  the  scale  of  horizontal  distance 
of  100  km.  This  indicates  that  the  advective  terms  are  not  negligible  in 
cases  when  the  horizontal  eddy  viscosity  terms  might  have  some  effects  on 
the  current  system. 


It.  TRANSPORT  OF  THE  CURRENT 

One  of  the  controversies  between  theories  and  observed  facts  of  the 
A.C.C,  is  the  relationships  between  the  zonal  transport  of  the  current 
and  the  wind  stresses.  According  to  the  theoretical  wd rks  of  Munk  and 
Palmln  (1951)  and  others,  due  to  the  wind  stresses  alone  the  total  zonal 
transport  becomes  larger  by  one  to  two  orders  of  magnitude  than  the  ob¬ 
served  value,  unless  large  frictional  forces  are  assumed  to  act  on  the 
current. 

In  order  to  estimate  the  zonal  transport  averaged  along  the  latitu¬ 
dinal  circles  surrounding  the  Antarctic  Continent,  the  averages  of  equa¬ 
tions  (1)  and  (5)  about  'A  from  0  to2/£are  taken.  Then  we  have 


Ur rRStnef  (U-e-  U.A  +  UAintrAuS/neVae 

+  RTa  J7/dJ  +  z£>-  <-osS  ir  =  -  ( 9  ft  si  n  &)'  (Pc  -  P,j 

__  _  (7) 

+  K  )ZZ  *  a  {R15i)'i9yl3(sines>u/30)/d9 

+  (zrcR-SnrP)  1  (9U/^/\)0~  /dA)^ 
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where  the  bar  indicates  the  average  along  a  latitude  circle.  The  suffixes 
0  and  1  indicate  the  values  respectively  at  the  east  and  west  coast  of  the 
Antarctic  Continent  or  South  America  where  the  latitude  circle  encounters 
it.  Therefore,  uG  -uq  and  p0  — pq  should  vanish  where  the  entire  latitude 
circle  lies  in  the  ocean.  When  the  coastlines  of  the  continents  are  con¬ 
sidered  parallel  to  meridians,  Uo  and  uq  vanish.  Although  the  coastlines 
may  not  consist  of  meridians,  u0  and  uq  are  approximately  equal  because 
the  coastlines  form  a  streamline  approximately.  Therefore,  the  terms 
including  u0  and  uq  in  equations  (7)  and  (8)  are  generally  negligible. 

If  the  current  along  the  coast  is  geostrophic,  p0  s  p-j_,  because  the 
coastlines  form  streamlines. 

The  values  of  (?>  j^and (3 are  not  equal  to  each  other.  How¬ 
ever,  the  order  of  magnitude  of  the  terms  including  the  boundary  values 
of  U  is 

A  AA  u6  (  ZTCF\*  Si  n  e  Ag)'1 


where  A  r  and  U  p  are  the  width  of  the  continent,  the  distance  of 

J  D  O 

the  coastal  boundary  layer  and  the  zonal  velocity  outside  the  coastal 
boundary  layer  respectively.  The  ratio  of  this  term  to  the  Coriolis 
term  may  become 

A  zxA  uB  (2/r  PAS'V©  ab-  fu)  ' 

This  equals  to  10“^  tfien  it  is  assumed  that  Q  A  30°) 


/o 


\ 3 


(21°  longitude).  Therefore,  the  terms  containing  the  boundary  values 
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of  0  l  ( jd/tnsy  be  negligible  compared  with  other  terms. 


The  integration  of  (7)  and  (8)  with  3  from  the  surface  to  the 

hr>+.+.nrn  vi  fll  rls  rp<;nof+.Ttrolir 


(11) 


averaged  along  a  latitude  circle  respectively,  H  is  the  depth  of  the  ocean 


and  C  ^  is  a  zonal  component  of  the  wind  stresses.  In  this  derivation, 


it  is  assumed  that  Uy  —  U(  vanishes  and  that  there  is  no  permanent  flux: 
of  water  across  the  sea  surface  or  the  bottom. 

Equation  (l)  can  be  integrated  immediately  and  the  solution  is 


(12) 


where  the  suffix  o  indicates  the  values  at  a  reference  co-latitude  taken 
probably  at  the  northern  limit  of  the  Antarctic  Circumpolar  Current, 

Munk  (1966)  estimates  that  the  bottom  current  of  the  effective  thickness 
1  km  (  s  h)  in  the  Antarctic  Ocean  has  northward  transport  with  the 
average  speed  of  0,07  cm/sec  (  Z  )•  The  Ekman  transport  in  the 

D 

upper  boundary  layer  is  southward  in  the  easterly  zone  near  the  Antarctic 
Continent  but  northward  in  the  westerly  zone  over  the  A.C,C,  Therefore, 
if  we  take  the  reference  latitude  at  the  boundary  between  the  east¬ 
erly  and  the  westerly  wind,  the  northward  transport  M0  consists  of  the 
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bottom  current  only  and  thus  equals  7  x  10^  cm2/sec.  When  the  water 
moving  northward  enters  the  three  oceans  bounded  by  the  continents,  the 
meridional  transport  h'i  may  be  modified  by  the  boundary  vaLues  of  [V! 
along  the  continents,  as  indicated  by  presence  of  the  u0  -u^  term  in 


A 


equation  (8).  If  the  meridional  transport  js‘/j  is  interpreted  at  the 
transport  below  the  upper  Ekman  layer  to  down  to  the  bottom  as  Munk 
(1966),  the  vertical  velocity  at  the  upper  boundary  does  not  necessarily 
vanish.  Then  the  northward  transport  from  the  Antarctic  Ocean  may  cause 
the  upwelling  in  the  mid-depth  of  the  central  part  of  the  three  oceans 
as  Munk  suggested* 

The  first  term  of  the  l.h.s,  of  equation  (10)  may  be  expressed  by 

(Vi  and  7.' 1  under  an  assumption  that  the  contribution  from  the  upper 

A  D 

Ekman  layer  to  'U^  is  neglected*  It  becomes 


r  H 


i 


o  £  Otsine)ck  - 


(13) 


The  second  term  of  the  l.h.s*  of  equation  (7)  may  be  negligible 
because  the  observed  data  indicate  d^//y£is  extremely  small  in  the  A.C.C. 
compared  with  the  one  in  other  oceans,  except  perhaps  in  the  surface  and 
bottom  boundary  layers  where  hj  almost  vanishes  and  thus  the  product 
also  becomes  small*  The  first  term  of  the  l.h,s*  of  equation 
(7)  can  be  determined  when  the  pressure  differences  at  both  sides  of 
the  continents  with  meridional  coastlines  is  determined.  The  barocline 
and  barotropic  pressure  difference  can  be  determined  by  hydrographic 
cast  and  by  measurement  of  the  mean  sea  level  differences  respectively. 
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However,  the  pressure  differences  may  be  small  if  the  geostrophic  equi¬ 
libriums  is  almost  valid  and  thus  this  term  is  also  dropped  in  the  fol¬ 
lowing  calculation* 

The  friction  at  the  bottom  can  be  determined  in  terms  of  /v  when 

f  A 

the  contribution  to  M  from  the  upper  Ekman  layer  is  neglected.  The 
horizontal  current  decreases  to  zero  at  the  bottom  within  a  small  depth 
of  the  bottom  Ekman  layer.  Therefore,  the  vertical  shear the 

bottom  can  be  represented  approximately  with  -  U,  O  j  where  ^and  D  are 

o  cr 

the  geostrophic  velocity  at  the  upper  boundary  of  and  the  thickness  of 

the  bottom  Ekman  layer  respectively.  The  geostrophic  current  Ug  is  ex- 

3  R  /  3  g 

•  Thus  the  friction  at  the  bottom  can  be 


H-O 


pressed  by  *  RPJ 

represented  by 

yK  s  u /qz I  ^  ~  K  i  Cco ^ ^  s  Pfess 


-l 


H -0 


H  . 

—  -  K  {9  R-f  D)  (3^/00) 


(ill) 


h-D 


where  y-  is  the  Coriolis  parameter.  The  thickness  of  the  bottom  Ekman 

I 

layer  D  equals  7l  (i  K. /£)■*•  Since  the  meridional  components  of  wind 
stresses  is  much  smaller  than  the  zonal  component  in  the  Antarctic  Ocean, 
the  contribution  of  the  Ekman  drift  to  may  be  negligible  compared 
with  the  geostrophic  component  of  [V]  as  seen  in  equation  (2).  Thus  we 
have 

p)cl X  —  f  '  He  k  (£  R)  j°/3 £  i 


/Ma-  [ 
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iSR)  (zp/du} 
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(15) 


where  He  is  the  equivalent  depth  if  the  pressure  gradient  at  the  whole 
column  is  proportional  to  its  value  at  the  depth  H-D.  When  the  pressure 
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gradient  is  barotropic,  He  is  exactly  equal  to  H.  Then  the  friction  at 
the  bottom  becomes 


jz  <tr 
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(16) 


where  ^  is  the  surface  elevation.  This  relation  is  the  same  as  obtained 
by  Kamenkovich  (1962)  except  the  numerical  constant  \)  z  7C  •  The  coef¬ 
ficient  of  in  his  expression  is  D/H  because  he  assumes  that 

D  1  •  When  the  relation  (l£)  is  substituted  in  (ill),  the 

friction  at  the  bottom  is  expressed  by 


K  ^  ^/o>  R 


K  (OHe)  lvlA=  -4  Ma 

(17) 


where  k  equals  about  3  x  10-7  (sec-T)  with  K  s  102  (cm2/sec),  He  z  3  Km 
and  D  =  10  m. 

By  use  of  such  expressions  of  different  terms  as  discussed  here,- 
equation  (11)  can  be  changed  into 
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0  S 1  o  0c 


H  R  S ,  n20  dQ 

=  ta- 
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The  second  term  of  the  l.h.s.  of  this  equation  equals  approximately 
-fh  }h  or  0.7  (cm2/sec2)  and  has  almost  the  same  orders  of  magnitude  as 
the  zonal  wind  stress'"/^  which  has  an  order  of  magnitude  of  1  cm2/sec2. 

v  A 

The  northward  transport  decreases  substantially  the  effect  of  the  stress 
of  the  westerly  wind.  The  total  eastward  transport  between  !?0°  and  60°  S 
is  about  l(P-h  crrP/sec  according  to  dynanic  calculation  of  Kort  (19^9) 
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based  on  reference  level  of  the  3000  decibar  surface*  Therefore.  V 

'  1 A 

10^  cm^/sec  and  thus  kM  ^  ~  0.3  (cm2/sec2).  Thus  the  bottom  friction  is 

also  of  the  same  orders  of  magnitude  as  , 

a  - 

If  it  is  assumed  that  the  Coriolis  term,  the  wind  stress  and  the 
bottom  friction  only  are  balanced  each  other  in  equation  (18),  we  have 


M 


X 


(Ta- 


n 

!  (A 


(19) 


where  T  represents  the  Coriolis  term.  Equation  (19)  yields  [A  “  IQlk 


A 


cm2/sec  for  T^  si  (cm2/sec2),  T0^S  0*7  (cm2/sec2)  and  k  :  3  x  10*"7 


(sec“A*  This  value  of  yields  the  eastward  transport  of  the  A.C.C* 
between  £0°  and  60°  S  equaling  10^  (cir^/sec)  which  is  exactly  of  the 
same  order  of  magnitude  as  observed  by  Kort  (195>9)*  Therefore,  the 
zonal  component  of  the  volume  transport  of  the  A.C.C,  depends  on  the 
meridional  transport  and  the  friction  coefficient  at  the  bottom,  as  well 
as  the  zonal  wind  stress  when  it  is  averaged  zonally.  The  horizontal 
eddy  viscosity  and  advective  terms  will  only  modify  the  zonal  distribution 
0fM>. 

In  order  to  determine  effects  of  the  horizontal  eddy  viscosity  and 
advective  terms  on  M  ,  equation  (18)  is  transformed  into 

A- 


( f?  +  h  i )  1  v-~  A 

(20) 


where  y  =  K|0is  the  northward  coordinate  and 

a,  =  A  cc- te/R  >  q,  =  Mc5,n^b  (Hs'h6)' 
-fe  (  =  p^lo5t^0o  Cot  Q  (  H  R  5 '  n  0) 

~j~\  -  2  -Cl  cc~t  & 


(21) 
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The  range  of  latitude  considered  is  between  6$°  and  35>°  S  or  2£°  ^  A0 
5>5>°.  The  values  of  MQ  was  based  on  estimation  by  Munk  (1966)  at  the 
north  of  the  Antarctic  convergence.  For  Mc  :  7  x  10®  (cm^/sec)  H  z  3  km 
A  =  10®  (cm^/sec)  and  0r  =  &5>°j  the  values  of  the  functions  of  (21)  at 
0  Z  23>°  and  5>J?°  are  listed  in  the  following  table. 

Table  1 

•  ai  (cm/sec)  a.2  (cm/sec)  kp  (sec-0  T7.  (cm^/sec®) 

0  =  25°  (65>°S)  34  x  10-1  2.0  x  lGr2  6.8  x  ICr^1  '  1.# 

0  =  $$o  (3^03)  1#i  x  1G-1  3.9  x  io-2  1**3  x  10-11  040 

This  table  indicates  that  a2  and  k^  are  negligible  compared  with  ap  and 
k  respectively,  suggesting  that  the  advective  term  with  the  meridional 
velocity  is  less  important  than  the  bottom  friction  and  horizontal  eddy 
viscosity  terms.  In  order  to  integrate  (21)  with  y,  a2  and  kp  are  ne¬ 
glected  and  ap  Is  expressed  by  a  term  inversely  proportional  to  y  for  a 
range  of  0)  from  2$°  to  £5>°.  If  equation  (21)  without  and  &p  Is 
represented  as  the  second  order  differential  equation  about  0  ,  the 
solution  can  be  expressed  by  Legendre  functions  of  complex  degree. 
However,  such  functions  are  inconvenient  both  for  analytical  and  numer¬ 
ical  treatment..  Therefore,  the  equation  will  be  reduced  to  one  type 
of  Bessel  equations  by  replacing  cot  0  in  ap  with  the  approximate 


expression.  Dimensionless  variables  UL  }>1  an<^  are  by  the 

relations :  ,  .  — r-  _ 

-  r'  A  =  lo  CV.* 

where  Mp,  L  and  T0  are  characteristic  scale  factors.  Equation  (20)  can 


be  transformed  into 


-t -wton&du/Jri-dZ  y.-nr\((y  (■>) 

'  ‘  (  t  (22) 
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where  n  =  i/R,  =  ft,  j_Z  ^  -? n  =  L *  To(M,A)’' 

The  function  cot  0  for  25°<10<:5&O  can  be  approximately  expressed  by 

L  ^  ^  with  the  appropriate  value  of  ('  and  range  of  /.  ,  which  may  be  deter¬ 
mined  so  as  to  yield  the  closest  agreement  between  the  two  functions. 

For  convenience,  the  approximate  function  is  set  to  be  exactly  equal  to 
cot  0  at  0  =  2^°  and  Q  s-  iiO°  and  the  value  of  ^  is  set  equal  to  one 

for  0  -  2f>°,  Then  ^  a  2,1 1|5  and  the  value  of  /?  corresponding  to  01  s 

!?5>0  becomes  2,6,  Two  curves  of  cot  and  Of1  are  plotted  in  Fig,  1, 
which  shows  the  difference  is  practically  insignificant  even  near  0  —  2,6, 
The  characteristic  distance  L  can  be  determined  as  2,1  x  10®  cm  from  the 
condition  that  the  co-latitudes  2f>°  and  3>5°  correspond  respectively  to 
/?  a  1  and  2,6,  Equation  (22)  is  rewritten  as 


M  o  a  f'/t'-  3e*-/C  -  on  (T-v-  70  (23) 

where  the  primes  indicate  differentiation  by  0  , 

A  forcing  function  7^ coming  from  the  Coriolis  term  is  proportional 

1  l/- 

too  cot  0  as  defined  in  (21),  and  its  dimensional  values  at  (j  a-  2^° 
and  f?5>°  are  listed  in  Table  1,  It  is  assumed  that  the  distribution  of 
another  forcing  function  ^  due  to  the  zonal  component  of  the  wind 
stresses  is  expressed  by 


to  Sin 


(2ii) 


When  the  westerly  belt  is  assumed  to  extend  from  65°  to  35°  S  '6=/T//.6 
/y  _  _  <  ,  The  characteristic  wind  stress  is  taken  as  1  cm^/sec^  and  the 

meridional  distributions  of  /^ancl  f  are  shown  in  Fig,  2  where  is  taken 

c  .  L  a  p 
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as  2  corresponding  to  the  maximum  wind  stress  of  2  cm^/sec^  at  the  latitude 
*0°  S. 

For  the  characteristic  values  of  A  s  10®  (cm2/sec),  L  a  2*1  x  10® 

(cm),  Mq_  ~  106  (cm^/sec),  /  a  (1  cm^/sec2)  and  k  -  3  x  lCr?  (sec~l), 
the  dimensionless  parameters  and  m  become  equal  to  132  and  Itlj.1 

respectively.  The  solution  of  equation  (23)  for  y>  ranging  from  1  to 

2.6  represents  the  meridional  distribution  of  M,^  from  latitudes  6£°  to 

2 

3?°  S.  Since  ^'and  m  are  much  larger  than  unity,  the  special  solution  of 
(2b)  can  be  given  approximately  by 


in  yd  2  ( 


(2S) 


This  represents  the  zonal  transport,  when  the  horizontal  eddy  viscosity 
terms  are  neglected.  It  should  be  noted  that  the  zonal  transport  becomes 
westward  from  65°  to  5>9°  S  in  case  of  'T' =  9  due  to  the  effect  of  the 

^  O 

meridional  transport  even  with  the  eastward  wind  stress,  as  seen  in 
Fig.  2.  The  elementary  solutions  of  a  homogeneous  equation  from  (23) 
are  expressed  by 
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respectively  the  first  and  second  kind  of  Bessel  function  of  the  order  Jk' 


/~~7  '  '  '  — ,  (  2. ) 

where  Cq_  and  C2  are  integration  constants  and  and  .  are 


equaling  (\~a~)/'y 

Because  yy  /  and  -2-6 *  the  Bessel  functions  (  K;ll) 

ry  1  2  ) 

and  L.  K  >1 1  )°an  ®e  expressed  with  their  asymptotic  formulas.  Owing 

\ 

to  the  imaginary  argument  ,  the  asymptotic  expansions  of  Hankel 
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functions  are  more  convenient  (Erdelyi  et  al,  195&)*  For  the  boundary 
condition  that 

yU.  —  O  oJr  J  4  q~  2.6  (27) 

the  solution  of  (23)  can  be  expressed  by  the  first  term  of  asymptotic 
expansions  as 

AL  -  /4  -  (?,  /'t ) "  % ,  1  -  x  A-  >b  )\  -  &  ^ 


(28) 

where  ^  and  are  respectively  the  higher  and  lower  latitude  boundaries 

(  n  -  \  n  -  2-  6  in  the  present  model)  and/-/  and//,  are  the  values 

L  t  ~  I  z  /  s  i  '  jZ. 

of y/ at  /^and  /^respectively.  The  exponential  terms  of  (28)  represent 
the  behaviors  of  /  ^  in  the  boundary  layers  near  IJ  i.  Ij  ^  and  ^  .  The 
thickness  of  the  boundary  layers  is  proportional  to  »  If  the  pro¬ 

portional  constant  is  taken  as  2.3,  the  thickness  of  the  horizontal 
boundary  layer  becomes  0.21  in  the  dimensionless  variable  //  or  1*20  km  in 
the  actual  distance  the  present  choice  of  values  A  and  kr.  Numerical 
calculation  of  the  solution  (28)  is  made  for  A  -  10®  (cm^/sec)  and  the 
distributions  of  Al  against  /?  and  ^0  are  plotted  in  Fig.  3.  The  result 
indicates  that  the  boundary  layer  solutions  at  -  /i  and  /?  do  not 
influence  the  transport  value  in  the  interior  of  the  Antarctic  Ocean. 
The  effect  of  the  spherical  shape  of  the  earth  is  represented  by  the 
coefficient  C^of  equation  (23)  and  the  solution  (28)  but  it  does  not 
change  the  feature  of  the  boundary  layers  very  much  because,}/ is  large 

and  thus  a  latitudinal  extension  the  boundary  layer  is  limited.  How- 
*  * 

ever,  if  ^  is  of  the  order  of  unity,  the  Hankel  functions  in  (26)  can 


no  more  be  represented  accurately  by  the  first  term  of  the  asymptotic 
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expansions.  This  indicates  that  the  thickness  of  the  boundary  layers  in¬ 
creases  to  the  order  of  a  few  thousand  kilometers,  and  thus  the  effect  of 
sphericity  of  the  earth  influences  in  the  entire  boundary  layer.  However, 
is  proportional  to  A“2  and  thus  A  must  be  at  least  of  the  order  of 
10-^  cm^/  sec  for  )£'  being  of  the  order  of  unity.  Although  such  large 
value  of  A  was  assumed  in  earlier  works  (Hidaka  and  Tsuchiya  for  instance), 
the  values  of  the  horizontal  viscosity  which  is  encountered  in  theories 
of  the  circulation  in  other  oceans  are  of  the  order  of  10®  cm^/sec. 

In  this  calculation  the  maximum  7^  of  the  wind  stress  function  is 

^  o 

taken  as  1*5  and  2,  corresponding  to  the  maximum  wind  stress  of  1.5  and 

2  cm2/sec2  respectively.  Fig.  3  indicates  that  the  maximum  transport  is 

3*8  x  10^  and  2.2  x  10^  (cm2/sec)  respectively  for  -  2  and  1.5*  The 

numerical  integration  of  the  transport  curves  of  Fig.  3  indicates  that 

the  total  eastward  transport  is  5*8  x  10^  (cm^/sec)  between  59*5°  and 

36.5°  S  and  2*8  x  10^  (cm2 /sec)  between  57*3°  and  38.3°  S  for  'T~  ~  2 

o 

and  1.5  respectively  and  that  the  transport  between  55°  and  1*5°  S  is 
3.7  x  10-^  and  1*9  x  10^  (cm^/sec)  for  '7'*  —  2  and  1.5  respectively. 
Therefore,  the  numerical  values  of  the  transport  for  ^  s  2  are  two  to 
three  times  larger  than  Kort*s  (1959)  value  calculated  from  the  ob¬ 
served  data.  The  important  point  is  that  the  meridional  transport 
drastically  decreases  the  eastward  transport  caused  by  the  eastward 
wind  stress  only.  (When  the  eddy  viscosity  A  is  taken  107  (cm  /sec), 
the  transport  also  decreases  almost  to  half  of  the  values  for  A  - 
10®  cm^ /sec  as  indicated  in  the  analysis  of  Appendix.  In  this  case  the 
special  solution  cannot  be  expressed  ty  equation  (25).  However,  choice 
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of  the  value  of  A  is  rather  arbitrary  except  for  its  orders  of  magnitudes* 
Therefore,  the  better  agreement  with  the  observed  data  does  not  neces¬ 
sarily  mean  that  the  larger  value  of  A  is  more  valid.) 

5,  DISCUSSION 

An  effect  of  the  meridional  transport  on  the  A.C.C,  seems  to  be 
very  important.  The  existence  of  the  general  northward  transport  is 
indicated  indirectly  from  conservation  of  salinity,  carbon- lli,  oxygen  and 
heat  in  the  deep  ocean  (Munk,  1966),  This  transport  may  be  caused  by 
sinking  of  the  upper  water  along  the  periphery  of  the  Antarctic  Conti¬ 
nent,  It  is  not  the  Sverdrup  transport  as  discussed  by  Wyrtki  (1962), 
because  it  does  not  depend  on  the  wind  stress  curl. 

In  order  to  determine  the  circulation  due  to  sinking  of  heavy 
water  along  the  Antarctic  Continent,  a  simple  experiment  was  made  with 
our  rotating  table,  A  circular  island  was  installed  at  the  center  of 
a  circular  tank  with  a  diameter  of  130  cm.  The  island  has  a  diameter 
half  of  that  of  the  tank  at  the  water  level  and  its  coast  sL  opes  towards 
the  periphery  of  the  tank  with  an  angle  of  !i3>°.  The  table  was  rotated 
with  constant  speeds  ranging  from  Ij.  rpmto  8  rpm.  Diluted  ink  with  cooled 
water  was  introduced  manually  on  the  periphery  of  the  island.  The  ink 
was  slightly  heavier  than  the  water  in  the  tank  and  flowed  down  the 
slope.  One  set  of  the  experiment  corresponding  to  the  rotation  speed 
of  £  r,p.m,  is  shown  in  Fig,  lu  The  rotation  of  the  tank  was  counter¬ 
clockwise,  corresponding  to  the  northern  hemisphere.  The  ink  spread 
clockwise  along  the  slope  and  at  the  bottom,  generating  an  eastward 
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It  should  be  remembered  that  the  zonal  and  meridional  transports 
discussed  here  are  those  averaged  along  latitude  circles.  When  a  lat¬ 
itude  circle  encounters  the  coastlines,  equation  of  continuity  becomes 


3 

a© 


sin  e)  +  HE 


(M 


AO 


M 


(29) 


instead  of  (11),  where  andf^^are  values  of  the  zonal  transport 
at  the  east  and  west  coast  of  the  continent  respectively.  If  the 
streamline  strictly  follows  the  coastline,  .  However,  the 

streamlines  in  our  model  experiments  indicate  a  weak  clockwise  eddy 
in  the  Weddel  Sea  and  a  wake  area  to  the  east  of  South  America  in 
agreement  with  the  transport  lines  obtained  by  Kort  (19^9)* 

Cne  example  of  the  streamlines  in  the  neighborhood  of  Drake  Passage 
determined  by  the  experiments  is  schematically  shown  in  Fig.  In  this 
experiment,  a  driving  force  was  the  wind  generated  by  eight  wind-blowers 
mounted  on  the  rotating  table  and  the  rotating  speed  of  the  table  was 
£  r.p.m.  In  such  circumstances  Nj  may  be  slightly  negative  (west- 

ward)  or  at  least  almost  vanishes  but  j^^is  positive  and  of  the  order 
of  10^  to  10^  cm^/sec.  Integration  of  equation  (29)  with  respect  to  0 
indicates  that £[n0.| 'J  may  increase  northwards  within  the  latitudes 
corresponding  to  the  Drake  Passage  by  about  U  x  10^  cm2/sec,  if  it  is 

assumed  that  M  s  10^  cm2/ sec.  This  increase  in  the  meridional  trans- 

'  'IX 

port  may  correspond  to  the  eastward  wind  stress  of  O.li  cm2/s  ec2.  Al¬ 


though  a  part  of  this  effect  is  cancelled  by  the  difference  of  the 
pressure  in  equation  (10 ),  the  effect  of  the  Drake  Passage  is  to 
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diminish  the  eastward  transport  through  increase  of  the  averaged  north¬ 
ward  transport*  More  quantative  evaluation  of  this  effect  will  be  made 
by  analysis  of  the  boundary  value  problem  and  by  experiments  with  the 
rotating  model  and  will  be  reported  later* 

In  the  later  experiments  the  periphery  of  the  model  Antarctic 
Continent  was  cooled  and  the  rim  of  the  tank  was  heated  to  simulate 
the  circulation  caused  by  differential  heating  in  the  meridional  di¬ 
rection*  The  details  of  these  experiments  will  be  described  in  the 
second  paper.  The  horizontal  circulation  in  these  experiments  shows 
three  clockwise  gyres  probably  owing  to  the  coastline,  instead  of  the 
continuous  circular  westward  flow  around  the  Continent  as  expected  from 
the  experiment  of  Fig*  li*  Also  in  these  experiments  the  meridional 
components  of  the  surface  and  bottom  current  were  generally  southward 
and  northward  respectively,  but  their  zonal  components  were  almost  in 
the  same  direction*  The  meridional  shearing  motion  may  be  due  to  the 
presence  of  the  outer  rim  in  the  experiments*  However,  predominance 
of  the  vertically  uniform  zonal  component  of  the  current  in  the  exper¬ 
iments  confirms  the  validity  of  assumption  of  the  bottom  friction 
proportional  to  the  zonal  transport* 
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EXPLANATION  OF  FIGURES 


Fig.  1  Comparison  of  cot  (3  (Qi  co-latitude,  Curve  A)  and  2.lli!?  r{~^ 
(Curve  B).  ^ 

Fig.  2  Meridional  distributions  of  T^the  zonal  component  of  the 

dimensionless  wind  stress 'Curve  A)  and  of  (the  dimensionless 
Coriolis  term,  Curve  B).  lr 

Fig.  3  Meridional  distributions  of  the  dimensionless  zonal  transport* 
(Curve  A  for  "T0  =  2,  Curve  B  for  T  =  1*5>  To  is  the  maximum 
value  of  the  dimensionless  zonal  wind  stress). 

Fig.  li  Schematical  pattern  of  dye  released  at  the  periphery  of  an 

island  with  a  sloping  coast  in  a  rotating  basin.  (The  rota¬ 
tion  rate  was  6  r.p.m.  The  pattern  was  drawn  based  on  16  ran 
motion  picture  film  seven  minutes  after  the  release  of  dye.) 

Fig.  5>  Schematical  presentation  of  the  streamlines  near  Drake  Passage 
based  on  the  rotating  model  experiment.  (The  rotation  rate  is 
5>  r.p.m.  Driving  forces  are  wind  only.) 

Fig,  A-l  Comparison  of  the  original  forcing  function  (Curve  A)  and  the 
approximate  one  (Curve  B). 

Fig.  A-2  A  meridional  distribution  of  the  dimensionless  zonal  transport 
for  A  a  IQ?  (cm2/sec). 
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APPENDIX 


A  particular  solution  of  the  non-homogenous  equation  (25>)  can  be 
obtained  as  a  polynomial  of  /7':if  the  forcing  function  is  expressed  by 

a  similar  polynomial.  For  instance,  when  7"'  —  ^  is  given  by 
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a  particular  solution  can  be  expressed  by 

/  A 


/^s  ~  ^  ^  z‘2  ^  ^4. 

Substituting  (A-l)  and  (A-2)  into  equation  (23)  and  equating  coefficients 
of  the  same  power  of  T)  to  zero,  we  have 
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Therefore,  tj~2n  an^  t  ?an  exPresse^  successively  in 

terms  of  Cir>  from  equations  (A-3)  (A— U )  and  (A-6),  An  exact 

particular  solution  like  (A-2)  is  more  suitable  for  a  case  of  A  z  10? 
(cm^/sec)  than  an  approximate  one  like  equation  (25)*  Therefore,  the 
function  7"  _  7-"  is  expressed  in  a  polynomial  up  to  H  ^  ,  The  constants 

y\  l1r  1 

C0,  C2  and  Cj^  are  determined  so  that  the  polynomial  may  have  the  same 
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values  as  the  original  function  of  >  -  r  at  south  and  north  boundaries 

A 

(  rj-  I  u-Atl  of  the  range  of  V  and  also  the  same  maximum  value  as 

the  original  function*  Comparison  between  the  polynomial  and  the  original 
function  is  shown  in  Fig*  A-l*  The  results  of  the  complete  solution  of 
equation  (2^)  satisfying  the  boundary  conditions  U-Q  at  O-  /  and  2*6 
are  shown  in  Fig*  A-2,  It  is  seen  that  the  value  of  A  3  10?  (cm^/sec) 
decrease  the  meridional  transport  almost  by  kOf°  from  the  one  for  A  z  10® 

(  cm^/sec)  as  stated  in  the  text.  The  total  eastward  transport  (with 
z  2)  becomes  3*6  x  10^  (cm?/sec)  compared  with  5>*8  x  10^  (cm®/sec)  for 
A  =  108  (cm^/sec).  The  transport  between  h%°  S  and  55°  S  becomes  2*0  x 
10-^  (cm3/sec)  compared  with  3.7  x  10^  (cm8/sec)  for  A  =  108  (cm^/sec). 
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